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SUMMARY 


This  document  contains  a  con-sis^pai  comprehensive  treatment  of 
the  current  knowledge  of  solar-particle  events. 

A  qualitative  description  has  been  given  of  the  various  methods  of 
detecting  solar-particle  events,  and  all  of  the  pertinent  data  available 
have  been  tabulated. 

The  characteristics  of  each  of  the  events  have  been  deteimined  and 
from  these,  the  intensities  and  the  characteristic  rigidities  have  been 
estimated  for  each  event. 

Based  upon  this  information,  the  radiation  doses  have  been  calcu¬ 
lated  for  each  event  and  cumulative  dose-probability  curves  have  been 
derived. 

Comments  are  included  regarding  missions  in  solar  cycle  20  and 
various  techniques  for  predicting  the  occurrences  of  solar-particle  events 
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INTRODUCTION 

Although  the  characteristics  of  the  galactic  cosmic  radiation  have 
been  recognized  and  investigated  in  detail  for  many  years,  the  ability  of 
the  sun  to  produce  large  bur^s  of  cosmic  radiation  has  been  studied  in 
detail  only  since  1956.  In  the  fifteen  years  preceding  the  great  solar- 
cosmic-ray  event  of  February  23,  1956,  there  were  only  four  cases  of  an 
increase  of  cosmic-ray  intensity  at  the  earth  (as  measured  at  ground  level 
by  ionization-type  detectors)  that  could  be  related  to  specific  solar 
activity  (Ref.  1).  Therefore,  the  idea  became  widespread  that  these 
events  were  quite  rare. 

In  the  period  subsequent  to  1956,  a  wide  variety  of  new  techniques 
have  been  employed  to  advance  our  understanding  of  the  solar-cosmic-ray 
events.  These  include  particle  counters  and  emulsions  flown  in  balloons, 
satellites,  and  space  probes,  as  well  cs  measurements  of  reduction  of 
cosmic  noise  coused  by  polar  cap  absorption  (PCA). 

With  these  techniques,  studies  of  solar-particle  bursts  over  the  last 

few  years  are  beginning  to  provide  a  reasonably  complete  picture  of  these 

events  (Refs.  2,  3,  and  4).  The  cosmic-ray  particles  ejected  from  the  sun 

are  known  to  be  primarily  protons  with  typically  steep  energy  (or  rigidity) 

spectra.  Particle  energies  range  from  less  than  10  Mev  to  a  few  Gev  on 

occasion.  Alpha  particles  and  heavier  nuclei,  principally  the  C,  N,  O 
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group,,  may  also  be  present  In  significant  numbers. 

These  solar  cosmic  rays  begin  to  arrive  near  the  earth  within 
minutes  after  they  are  produced  in  the  region  of  a  flare  on  the  sun.  The 
intensity  of  the  solar  cosmic  rays  is  a  strong  function  of  solar  activity.  Such 
cosmic-ray  particles  are  present  near  the  earth  in  detectable  numbers  a 
significant  fraction  of  the  time  (5  to  10  percent)  and  exceed  the  normal 
galactic-cosmic-ray  intensity  fora  much  smaller  period  (~1  percent). 
Occasionally,  for  intervals  of  one  or  two  days,  their  intensities  are  many 
thousands  of  times  greater  than  the  galactic-cosmic-ray  intensity.  On  the 
other  hand,  at  periods  near  sunspot  minimum,  one  to  two  years  may  elapse 
between  significant  solar-particle  outbursts. 

The  intensity-time,  directional,  spectral,  and  charge  characteristics 
of  solar  cosmic  rays  differ  considerably  from  event  to  event.  For  example, 
most  events  have  many  more  low  energy  particles  than  high  energy  ones. 

But  some  events,  which  produce  a  rare  effect  in  a  sea-level  detector, 
contain  a  relatively  large  fraction  of  high-energy  particles,  correspondingly 
few  low-energy  particles,  and  thus  give  a  small  integrated  particle  intensity. 
This  latter  characteristic  causes  the  event  to  be  classified  as  a  relatively 
small  one.  Other  important  dissimilarities  exist  from  event  to  event  and, 
for  this  reason,  we  shall  catalog  some  of  the  more  important  ones  individually. 

In  contrast  to  the  solar  cosmic  rays,  the  galactic  cosmic  radiation  , 
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exhibits  a  rather  regular  behavior,  it,  too,  consists  primarily  of  protons 
with  a  significant  admixture  of  alpha  particles  and  heavier  nuclei.  The 
proportions  of  the  various  charge  components  remain  constant  throughout 
the  cycle  of  solar  activity  while  the  total  particle  intensity  changes  by  a 
factor  of  five.  Maximum  intensity  of  galactic  cosmic  rays  occurs  at 
sunspot  minimum  in  anti-correlation  with  the  solar  activity.  The  spectrum 
changes  only  slowly  during  the  course  of  the  sunspot  cycle.  Near  sunspot 
minimum  there  are  mora  low-energy  particles  and,  consequently,  the 
spectrum  is  steeper.  In  fact,  primarily  the  lower-energy  particles  (<  5  Gev) 
are  affected  as  c  result  of  the  solar  activity. 
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CHARACTERISTICS  OF  SOLAP-COSMIC-RAY  EVENTS 

Before  considering  some  of  the  feotures  of  the  production  of  solar 
cosmic  rays,  we  shall  establish  some  sort  of  workable  picture  of  the 
characteristics  of  the  solar-cosmic-ray  events  themselves.  As  we  have 
pointed  out,  the  details  of  each  cosmic-ray  event  may  differ  considerably. 
Nevertheless,  we  may  define  certain  general  characteristics  common  to 
all  events. 

INTENSITY-TIME  PROFILE 

A  typical  event  will  hove  a  set  of  intensity-time  profiles,  one  for 
each  energy,  such  as  those  shown  in  Fig.  1,  Certain  times  characterize 
each  profile. 

Onset- D elay  Time 

The  onset-delay  time  is  defined  as  the  time  from  the  maximum  intensity 
of  the  optical  flare  to  the  arriW  of  the  first  particles  at  the  earth.  This 
time  is  variable  from  event  to  event  and  is  strongly  energy-dependent,  with 
the  higher  energies  arriving  first.  It  also  depends  on  the  location  of  the  fiore 
on  the  sun,  with  those  on  the  west  portion  of  the  visible  sun  having  shorter 
onset-delay  times  than  those  on  the  east  portion.  The  onset-delay  times  may 
vary  from  a  few  minutes  for  high-energy  particles  (>  500  Mev)  in  some  events 
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to  many  hours  for  low-energy  particles  (<  30  Mev)  in  other  events. 

Rise  Time 

The  rise  time  is  defined  as  the  time  interval  at  the  earth  from  the 
first  arrival  of  particles  of  a  particular  energy  to  the  time  at  which  the 
maximum  intensity  of  these  particles  is  attained.  This  time  is  strongly 
energy-dependent,  for  the  higher  energies  reach  maximum  intensity  first. 
There  is  evidence  that  this  rise  is  exponential  with  time  (Ref.  5)  in  a  num¬ 
ber  of  events.  This  allows  us  to  define  a  characteristic  rise  time,  t^.  This 
time  is  usually  closely  related  to  the  onset-delay  time  in  a  particular 
event  and  depends  in  a  similar  manner  on  the  location  of  the  flare  on  the 
sun.  The  rise  time  may  range  from  a  few  minutes  for  high-energy  particles 
in  some  events  occurring  near  the  west  limb  of  the  sun  to  many  hours  for 
low-energy  particles  from  events  occurring  near  the  east  limb  of  the  sun. 

Decay  Time 

A  growing  body  of  evidence  (Reis.  3  and  5)  indicates  that  the  decay 
of  particles  above  a  particular  energy  (or  rigidity)  is  also  exponential  at 
most  rimes  and  over  a  wide  range  of  energies.  Thus  we  may  define  a 
characteristic  decay  time,  tp,  and  write  the  following  equation  for  the 
intensity  of  flare  particles  with  energies  greater  than  some  energy,  E,  at 
some  time,  t,  after  the  maximum  intensity  of  particles,  lmax(IE),  is  attained: 
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The  time  t^,  a  function  of  energy,  is  smaller  for  higher  energies.  It  is 
variable  from  event  to  event,  although  the  variation  seems  to  be  influenced 
more  by  general  interplanetary  conditions  than  by  the  location  of  the  flare, 
which  governs  the  onset-delay  and  rise  times.  A  sudden  change  in  the 
interplanetary  magnetic  field  configuration  during  the  course  of  an  event 
may  drastically  affect  the  decay  rate,  as  measured  by  the  decay  time.  In 
simple  events,  this  time  may  range  from  two  to  four  hours  for  high-energy 
particles  to  two  to  four  days  for  low-energy  particles.  Generally,  then, 
the  decay  time  is  longer  than  the  rise  time,  in  some  instances  up  to  ten  times 
as  long. 

ANISOTROPIES  IN  SOLAR  COSMIC  RAYS 

For  the  most  part,  the  solar  flare  particles  reach  the  vicinity  of  the 
earth  very  nearly  isotropically  (within  ±  5  to  10  percent),  and  the  afore¬ 
mentioned  ir.tensity-time  characteristics  refer  to  this  isotropic  radiation. 
However,  significant  anisotropies  in  the  direction  of  arrival  of  solar-flare 
particle;  at  the  earth  do  exist.  In  particular,  the  particles  appear  to  arrive 
from  a  highly  preferred  and  fairly  narrow  direction  30°  to  60°  west  of  the 
earth-sun  line  (Ref.  6).  This  behavior  is  most  pronounced  for  the  high- 
energy  particles.  This  directional  radiation  is  usually  of  short  duration 
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compared  to  the  isotropic  part,  but  intensities  about  10  to  100  times  the 
maximum  isotropic  radiation  may  occur  in  the  preferred  direction  fora  few 
minutes.  As  a  result  of  the  action  of  the  geomagnetic  field,  this  radiation 
approaches  the  surface  of  the  earth  in  certain  well-defined  areas  called 
impact  zones. 

ENERGY  SPECTRUM  OF  FLARE  PARTICLES 

It  is  obvious  from  the  foregoing  discussion  of  the  intensity-time 

characteristics  of  the  flare  particles  and  from  the  energy  dependence  of 

each  of  the  parameters  involved  that  the  measured  energy  spectrum  for  any 

one  event  will  depend  upon  the  specific  time  during  the  event  at  which 

the  measurements  are  made.  The  problem  of  determining  o  useful  and 

meaningful  representation  of  the  solar  flare  particle  spectrum  has  been 

one  of  the  most  difficult  connected  with  the  study  of  the  flare  particles. 

Usually,  the  differential  number  spectrum  of  the  flare  particles  has  been 

expressed  either  as  a  power-law  energy  spectrum,  dJ/dE  =  K^/E  or  as 

n2 

a  power-law  rigidity  spectrum,  dJ/dP  -  K^/lP  (the  rigidity  of  a 

particle  -  P  =  pc/2e,  where  p  is  the  particle  momentum,  c  the  velocity 

of  light,  e  the  unit  charge,  and  Z  the  charge  number  of  the  particle). 

When  either  of  these  representations  is  used,  n^  (or  n£)  is  a  function  of 

both  energy  (or  rigidity)  and  time.  For  example,  n^  (or  ng)  is  usually 

smaller  for  lower  energies  or  rigidities  because  of  the  "bending  over"  of 
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the  spectrum.  Also,  (or  n increoses  with  time  after  the  onset  of  an 
event,  as  the  particle  spectrum  steepens  with  time. 

Because  the  energy  or  rigidity  interval  viewed  in  any  one 
measurement  is  usually  quite  narrow,  it  is  possible  to  represent  reasonably 
well  the  number  spectrum  of  solar-flare  particles  as  a  power-law  energy 
or  rigidity  spectrum  with  an  n  which  is  constant  over  the  interval  of 
measurement.  This  does  not  prove  to  be  a  useful  representation  over  the 
entire  range  of  rigidities  involved  in  the  solar  outburst,  however. 

After  a  careful  examination  of  the  data  from  a  large  number  of 
events,  which  were  difficult  to  interpret  if  the  solar-particle  spectrum 
were  considered  to  be  a  power  law  in  either  energy  or  rigidity,  the  con¬ 
clusion  has  been  reached  that  an  exponential  rigidity  spectrum  best 
describes  the  solar  flare  particles  (Ref.  4).  This  spectrum  is  of  the  form 
dJ/dP  =  dJ/dP^  expt-P/P^t)].  The  characteristic  rigidity  Pq  in  a 
specific  event  is  a  function  of  time  only,  and  Pq  decreases  as  the  event 
progresses.  Its  value  and  time-dependence  vary  from  event  to  event. 

Such  spectra  apply  to  all  energies  greater  than  2  Mev  if  the  time  from  the 
flare  is  a  few  times  the  onset-delay  time  for  particles  of  a  particular 
energy. 

The  time  development  of  a  typical  solar-cosmic-ray  spectrum 

during  an  event  is  shown  in  Fig.  2.  Typical  values  of  the  characteristic 
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FIGURE  2.  Typical  energy  spectra  obtained  at 
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rigidity  P  .  associated  with  times  near  the  occurrence  of  I  for  the 
low-energy  particles,  may  range  from  about  300  Mv  in  events  with  very 
flat  spectra  containing  many  high-energy  particles  to  about  50  Mv  in  events 
with  very  steep  spectra  containing  many  low-energy  partic'es  and  producing 
strong  radio-wave  absorption. 

CHARGE  COMPOSITION  OF  FLARE  PARTICLES 

Until  recently,  it  was  thought  that  the  particles  emitted  in  a  solar 

flare  burst  were  predominantly  protons.  Recent  observations  have  established 

the  presence  of  alpha  particles  in  at  least  nine  events  and  heavier  particles 

(principally  the  M  group  consisting  of  carbon,  oxygen,  nitrogen  and 

fluorine)  in  three  events  (Refs.  7  and  8).  These  heavier  particles  appear  to 

have  exponenrial  rigidity  spectra  similar  to  the  proton  spectra.  Thus,  at 

any  time  the  same  P^  applies  to  protons,  alphas,  and  heavier  nuclei.  The 

ratios  of  the  values  of  Jq  for  each  charge  group  (e.  g.,  J^/J  ,  J 

are  approximately  constant  during  nn  event,  although  this  ratio  varies  from 

event  to  event.  Observations  at  rigidities  greater  than  0.5  Gv  have  revealed 

a  J0p/-^0a  fhat  varies  from  one  to  about  twenty  with  a  predominance  of 

values  near  one.  The  ratio  J^^/J  deduced  from  observations  in  the 

same  rigidity  range,  appears  to  be  a  constant  from  event  to  event  with  a 

value  approximately  60.  Since  all  charge  components  appear  to  have  similar 

rigidity  spectra  above  ~  0.  5  Gv,  and  if  this  similarity  continues  to  low 
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rigidities,  then  in  events  where  J  /J  Is  about  one.  we  should  exnect 

op  on 

equal  totai  numbers  of  pi o tor.;  and  alpha  particles.  The  presence  of  these 

alpha  particles  has  not  been  directly  determined,  but  their  presence  may 

be  inferred  from  radio-absorption  measurements  (Ref.  9). 

If  an  alpha  particle  and  o  proton  have  the  same  range,  then  the 

alpha  particle  always  has  a  greater  rigidity,  exactly  double  at 

non-relativiitic  energies.  (Note  that  alphas  and  protons  of  the  same  range 

have  the  same  energy  per  nucleon. )  On  this  basis,  we  can  account  for 

the  experimental  observation  of  fewer  alphas  than  the  J  /J^  rot'0  pre- 

diets.  The  ratio  of  alphas  to  protons  with  a  given  range  increases  with 

decreasing  range.  Also,  this  ratio  is  a  function  of  the  characteristic  P^ 

2 

so  that  at  ranges  less  than  1  gm/cm  and  for  flat  spectra  (P^  >  100  Mv), 

the  ionization  produced  by  clphas  may  dominate  in  the  cases  where 

J  /J  is  about  one.  The  same  considerations  apply  to  the  heavier 
op  or  rr  ' 

nuclei  (principally  oxygen)  so  that  the  ionization  produced  by  these  par¬ 
ticles  can  be  a  significant  fraction  (~50  percent)  of  that  produced  by 

2 

protons  at  depths  less  than  1  gm/cm  . 


20 


D2-90469 


SOLAR  FEATURES  ASSOCIATED  WIIH  uOSMiC-RAY  OUTBURSTS 

If  is  clearly  established  that  solar  cosmic  rays  detected  near  the 

earth  originate  in  connection  with  large  flares  on  the  visible  disk  of  the  sun. 

Flares  are  sudden  increases  in  intensity  of  optical  radiation  over  large  areas 

in  the  solar  atmosphere.  The  size  scale  for  flares  is  based  on  their  area, 

with  the  largest  flares  (3+)  having  an  area  greater  than  1000  x  10  ^  of  the 

solar  surface.  Increases  in  the  optical  radiation,  which  occur  principally 

in  hydrogen  line  emission,  e.  g. ,  H  ,  form  a  negligible  fraction  of  the 

a 

over-all  energy  output  of  the  sun.  in  addition,  large  increases  in  solar 

energy  output  occur  at  both  ends  of  the  electromagnetic  spectrum:  at  wave- 
0 

lengths  less  10  A,  corresponding  to  X  rays  and  gamma  rays,  and  at  the  radio 
end  of  the  spectrum.  Of  particular  interest  in  association  with  the  cosmic- 
ray  output  are  the  radio  emissions  in  the  so-called  microwave  bursts  at 
frequencies  >3000  h/.c/s  (wavelengths  <  10  cm).  These  will  be  discussed 
later  in  greater  detail. 

A  flare  is  not  a  separate  or  distinct  event  but  originates  in  so-called 
active  regions  or  centers  on  the  sun.  Therefore,  ihe  study  of  flares  starts, 
in  effect,  with  observations  related  to  the  life  history  of  an  active  region. 
An  active  region  first  puts  in  its  appearance  as  a  small  plage  region  located 
in  one  of  the  two  sunspot  belts,  either  north  or  south  at  a  latitude  of  5°  to 

40°.  A  few  days  later,  sunspots  appear  in  the  region  and  proceed  to  grow 

21 


D2-90469 


in  number  and  in  size  as  the  size  and  brightness  of  the  plage  increases.  Sunspots 
ore  cooler  than  the  surrounding  photosphere  and  have  strong  magnetic  fields 
associated  with  them.  The  darker  central  region  of  a  spot  is  called  the  umbra, 
and  the  gray  region  around  it,  the  penumbra.  As  the  spot  group  continues  to 
grow,  small  flare s  occur  end  active  bright  filaments  (or  prominences)  begin  to 
appear.  When  the  spot  group  reaches  maximum  size,  the  flares  become  larger 
and  more  numerous.  At  this  stage  of  development,  or  somewhat  later,  cosmic- 
ray  events  are  most  likely  to  occur.  This  stage  may  be  reached  in  three  to 
four  days  in  some  instances,  although  typically  the  time  is  about  thirty  days. 

The  features  associated  with  the  cosmic-ray  flares  and  their  centers 
are  very  complex  indeed;  however,  certain  characteristics  appear  to  stand  out. 
These  are  described  below. 

CONNECTION  WITH  LARGE  FLARES 

Cosmic-ray  events  generally  occur  in  connection  with  the  larger  solar 
flares.  These  flares  have  areas  of  bright  emission  which  are  comparable  to  the 
area  of  the  entire  associated  sunspot  group  (~  2000  x  10^  of  the  solar  hemis¬ 
phere). 

LARGE  AND  COMPLEX  SPOT  GROUPS 

The  spot  group  itself  is  generally  large  and  complex  with  a  complex 
magnetic  field  configuration.  Also,  a  high  level  of  flare  activity  and 
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frequent  radio  bursts  are  associated  with  the  spot  group. 

H  RADIATION 
a 

Those  flares  that  produce  cosmic  rays  are  characterized  by  an 

unusual  flash  of  rudiation  in  H  .  It  Is  believed  that  during  this  period  of 

a 

one  to  two  minutes  the  cosmic  rays  are  accelerated  in  the  region  of  bright 
filaments  that  form  over  sunspot  umbrae.  Generally,  there  are  two  bright 
parallel  filaments  separated  by  the  magnetic  axis  which  Is  a  line  dividing 
the  main  North  poles  from  the  main  South  poles  within  the  group  (Ref.  10). 

MICROWAVE  RADIO  EMISSION 

The  flares  that  produce  cosmic  rays  are  also  characterized  by  micro- 
wave  (cm)  radio  emission  of  exceptional  intensity.  This  emission  appears 
and  disappears  with  a  very  rapid  time  scale  of  the  same  order  as  the  optical 
emission.  The  source  location  is  closely  associated  with  the  flare  region; 
the  height  of  the  source,  about  0.05  solar  radii  above  the  photosphere,  is 
constant  in  time  (Ref.  11).  The  source  has  roughly  the  same  dimensions  as 
the  flare  region.  Undoubtedly,  this  emission  ‘.s  the  radio  analogue  of  the 
optical  flash  and  is  an  essential  accompaniment  of  the  generation  by  these 
flares  of  relativistic  particles.  The  close  association  between  this  portion 
of  the  radio  emission  and  synchrotron  emission  of  electrons  moving  in  strong 
local  fields  is  an  Interesting  measure  of  the  particle  acceleration  (Ref.  12), 
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METHODS  OF  DETECIiNG  SOLAR  COSMIC  RAYS 

In  order  to  interpret  solar-cosmic-ray  events,  we  must  understand, 
among  other  things,  the  methods  used  to  detect  these  particles.  Before 
discussing  the  cosmic-ray  events  themselves,  then,  we  shall  briefly  discuss 
some  of  the  detectors  and  their  relative  sensitivities. 

NEUTRON  MONITOR 

The  neutron  monitor  consists  of  a  structure  of  BF^  proportional 
counters  surrounded  by  an  arrangement  of  lead  and  paraffin.  A  high- 
energy  (>  100  Mev)  neutron  or  proton  passing  through  the  lead  has  a  high 
probability  of  undergoing  an  interaction  in  which  some  of  the  reaction  pro¬ 
ducts  are  neutrons.  These  neutrons  are  slowed  down  in  the  paraffin  by 
collisions  with  hydrogen  nuclei  and  are  then  detected  by  the  proportional 
counters.  Since  neutron  monitors  are  usually  located  at  or  near  sec  level, 
the  detected  particles  are,  for  the  most  part,  the  secondary  nucleonic 
component  of  the  primary  radiation.  The  relation  between  this  secondary 
radiation  and  the  particles  incident  on  the  top  of  the  armosphere  is  a 
complicated  function.  It  is  the  result  of  a  nucleonic  cascade  and  is  best 
expressed  in  terms  of  a  specific  yield  function.  An  example  of  the 
relative  yield  of  secondaries  at  sea  level  os  a  function  of  rigidity  of  the 

primary  particles  at  the  top  is  shown  (Ref.  13)  in  Fig,  3.  It  can  be  seen  from 
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ii'iii  figure  that  the  neutron  Jenifer  is  relotively  insensitive  to  the 
less-fhon-1-Gv-rigidity  particles  which  form  the  main  bulk  in  events 
with  steep  spectra.  Only  events  where  substantial  numbers  of  particles 
above  1  Gv  occur  (P  generally  greater  than  100  Mv)  will  be  observed  by 
these  detectors  operating  near  sea  level. 

DIRECT  PRIMARY  PARTICLE  DETECTORS 

The  sensitivity  to  low-energy  particles  for  this  group  of  detectors 
Is  generally  limited  by  two  things:  (1)  the  geomagnetic  latitude  at  which 
the  observation  is  made,  and  (2)  the  amount  of  air  mass  above  the  detector. 
The  limitations  on  energy  as  a  function  of  air  mass  and  latitude  are  shown 
in  Fig.  4.  Since  the  essential  form  of  all  of  these  detectors  can  be  carried 
on  either  a  balloon,  rocket,  or  satellite,  the  difference  in  sensitivity  is 
mainly  one  of  air  mass  or  absorber.  Here  a  rocket  or  satellite  has  a 
distinct  advantage  over  a  balloon.  The  balloon-borne  detector,  of  course, 
has  an  overwhelming  advantage  over  a  neutron  monitor  at  sea  level  in  the 
detection  of  low-energy  particles. 

The  types  of  detectors  classified  as  direct  particle  detectors  include 
the  ones  described  below. 

Geiger  Counter 

This  device  simply  counts  the  total  (integral)  number  of  particles 
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passing  through  it.  Since  it  always  must  observe  solar  cosmic  rays  rolative 
to  the  background  of  galactic  cosmic  rays,  this  background  sets  the  limit 
to  the  sensitivity.  The  geiger  counter  may  be  used  to  give  time  and  direction 
discrimination  and,  with  absorbers,  crude  energy  measurements. 

Ion  Chamber 

This  device  counts  the  total  (integral)  ionization  of  the  particles  passing 
through  it.  It  also  must  observe  the  solar  cosmic  rays  relative  to  the  galactic- 
cosmlc-ray  background  which  sets  the  limit  to  its  sensitivity.  However,  since 
the  specific  ionization  of  the  solar  particles  is  greater  than  that  of  the 
galactic  ones,  an  ion  chamber  has  a  greater  sensitivity  than  a  comparable 
geiger  counter.  The  ion  chamber  may  be  used  to  give  time  measurements, 
but  direction  and  energy  measurements  using  it  alone  have  not  been  tried. 

It  has  been  used  to  advantage  in  conjunction  with  a  geiger  counter  to  obtain 
energy-spectrum  and  charge-composition  measurements  of  the  solar  particles. 

Energy-Loss  Measuring  Devices 

These  devices,  used  individually,  provide  information  similar  to  that 

obtained  by  using  an  ion  chamber,  but  they  may  be  readily  combined  in 

various  ways  for  use  in  telescopes  which  measure  directly  the  energy  and 

cnarge  of  the  particles  passing  through.  Such  detectors  include  scintillation 

2 

and  Cerenkov  counters  for  acceptance-solid-angles  (ASA)  exceeding  1  ern 

28 


D 2-90469 


sterodian  and  solid  state  detectors  for  ASA's  less  than  this  value.  Such 
devices  also  provide  time  and  direction  resolution.  Their  sensitivity 
exceeds  that  of  an  ion  chamber  and  geiger  Counter  under  comparable 
conditions  since  energy-loss  measuring  devices  can  observe  the  differential 
energy  spectrum  in  the  energy  region  of  response.  In  this  region,  the 
solar  cosmic  rays  can  often  be  distinguished  from  the  galactic  by  their 
rapid  decrease  in  number  with  increasing  energy.  Also  in  this  category 
are  nuclear  emulsions  which  provide  energy,  charge,  and  direction 
resolution  and  sensitivity  equal  to  or  surpassing  that  of  the  counters. 

Emulsions  do  not  provide  the  time  resolution,  which  Is  important  in  many 
aspects  of  solar-flare-particle  study.  Also,  the  need  to  recover  emulsions 
limits  their  use  to  balloons  and  local  rockets. 

INDIRECT  DETECTORS 

Most  indirect  detectors  respond  to  the  ionization  and  resul'ant 

radiowave  absorption  produced  by  the  solar  cosnrrc  rays  in  the  40-  to  100-km 

level  in  the  ionisphere.  Although  at  this  level  the  equilibrium  electron 

density  resulting  from  cosmic  rays  is  smaller  than  in  the  E  and  F  layers,  the 

large  collision  frequency  makes  the  product  of  collision  frequency  and 

density  sufficient  to  produce  appreciable  absorption  even  at  frequencies  as 

high  as  30  Mc/s.  Some  of  the  techniques  used  to  measure  the  effects  of 

this  ionization  are  described  below. 
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Measurement  of  Cosmic  Noise  Absorpt i on 

A  self-baiancing  receiver  called  a  rlomeiei  is  used  to  measure  the 
signal  strength  of  extraterrestrial  radio  noise.  The  increased  electron 
density  caused  by  the  solar  cosmic  rays  absorbs  the  extraterrestrial  noise 
more  effectively  and  consequently  decreases  the  signal  measured  at  the 
ground. 

In  Fig.  5  (Ref.  4}  the  total  absorption  at  30  Mc/s  is  shown  as  a 
function  of  the  spectrum  end  intensity  of  the  solar  particles  (assuming  pro¬ 
tons  only)  for  values  of  the  effective  cutoff  rigidity  appropriate  to  locations 
of  various  observing  stations.  Using  the  absorption  measurements  from  a 
single  station,  we  cannot  determine  either  the  spectrum  or  the  intensity 
of  the  solar  particles.  If  we  have  other  intensity  measurements,  e.g., 
direct  particle  measurements,  then  we  may  determine  both  the  spectrum 
and  the  intensity  as  precisely  as  we  can  calculate  the  absorption  due  to  a 
spectrum  of  solar  particles.  The  main  uncertainty  here  lies  in  the 
ionospheric  parameters  and  processes  that  determine  the  equilibrium 
electron  density.  Recent  advances  in  our  understanding  of  these  parameters 
and  processes  have  reduced  the  uncertainty  in  intensities  calculated  in 
this  manner  from  a  factor  of  about  ten  to  a  factor  of  two  or  less  (Ref.  14). 

The  spectrum  of  solar  particles  can  be  determined  by  use  of 
riometers  located  at  several  different  latitudes  as  well  as  by  use  of 
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riometers  operating  at  several  different  ficqu  cedes  r*t  a  single  location. 

There  is  a  difficulty  in  the  interpretation  of  riometer  data  because 
ionixation  produced  by  particles  not  directly  related  to  the  solor  cosmic 
rays  also  absorbs  extraterrestrial  radio  noise.  The  absorption  responsible 
for  this  problem  is  created  mainly  by  auroral  particles.  This  auroral 
absorption  is  generally  less  regular  in  both  its  temporal  and  spatial  charac¬ 
teristics  and  hence  is  separable  from  solar-cosmic-ray  absorption. 

Measurement  of  f  . 

_ _  _ mm 

In  this  technique,  a  radio  station  determines  the  minimum  signal 

frequency  which  can  be  sent  out  from  the  earth,  reflected  at  vertical 

incidence,  and  received  back  at  the  sending  station.  This  signal  is 

generally  reflected  from  the  E  or  F  layer.  Thus,  a  measurement  of  f  . 

is  essentially  a  measure  of  the  absorption  taking  place  at  the  lower  levels 

and  is  closely  related  to  the  results  of  absorption  measurements  by  riometers. 

A  large  network  of  stations  make  f  .  measurements.  In  some  instances, 

min 

the  observations  extend  back  to  1938.  Hence,  the  record  of  these 
observations  provides  one  of  the  best  available  monitors  of  the  occurrence 
frequency  of  solar-cosmic-ray  events,  A  quantitative  basis  of  the  f  . 
determinations  is  difficult  to  find  since  each  station  may  hove  o  somewhat 
different  signal  detection  limit.  Nevertheless,  we  show  in  Fig.  6  a  typical 
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relotion  between  the  absorption  of  30  Mc/s  and  the  change  in  f  .  from 

r  min 

its  normal  value  in  the  quiet  ionosphere.  Normal  values  of  f  .  are 

mm 

usually  less  than  1  Mc/s.  v; 

The  attenuation  of  the  f  .  signal  due  to  Ionization  below  the 
min 

reflecting  layer  is  plotted  as  a  parameter.  Curves  for  exponential  solar- 

particle  rigidity  spectra  with  average  rigidities  =  50  and  100  Mv  are 

shown.  It  is  apparent  that  Af  .  is  not  a  strong  function  of  P  .  Experi- 

min  o 

mental  measurements  of  Af  .  and  of  the  absorption  at  30  Mc/s  In  one 

mm  r 

event  are  shown  by  the  open  circles  in  Fig.  6  (Ref.  15).  These  lie  along 

a  curve  which  suggests  for  the  f  .  receiver  a  detection  capability  of 

min 

-60  db  in  this  instance.  Furthermore,  this  curve  has  an  average  slope  of 

abour  1  Mc/s  change  in  f  .  for  each  0.6  db  absorption  at  30  Mc/s. 

mm 

Hence,  this  method  of  measuring  is  comparable  in  sensitivity  to  riometer 

measurements  of  cosmic  noise  absorption.  Measurements  of  f  .  could  be 

min 

improved,  however,  by  increasing  the  sensitivity  of  the  receivers  or  by 
actuclly  measuring  the  strength  of  the  returning  signal. 

By  this  method,  just  as  by  the  riometer  methods,  ionization  due 
to  auroral  particles  can  be  mistaken  for  that  produced  by  solar-cosmic-ray 
events. 


Forward-Scatter  Techniques 

In  this  method,  pioneered  by  Bailey  (Ref.  16),  a  signal  is  reflected 
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at  grazing  incidence  from  the  E  or  F  layer.  Reflection  from  a  shallow 

angle  of  incidence  is  equivalent  to  scattering  through  a  small  angle, 

hence  the  term  forward  scattering.  The  attenuation  of  the  remotely 

received  signal  is  a  measure  of  *he  absorption  in  the  reflecting  layer. 

Tii is  method  is  more  sensitive  to  lav/  levels  of  absorption  because  of  the 

very  long  signal  path  through  the  absorbing  layer.  Also,  it  effectively 

discriminates  against  spatially  localized  auroral  absorption.  However, 

this  method  has  not  as  yet  been  put  on  a  quantitative  basis. 

In  conclusion,  in  Table  1  we  compare  the  capabilities  of  these 

various  methods  of  detection  for  spectra  with  P  =50  and  100  Mv.  For 

o 

each  method,  we  assume  the  smallest  detectable  response  due  to  solar 

cosmic  rays  and  translate  this  into  a  corresponding  flux,  J  ,  of  solar 

o 

cosmic  rays.  The  integral  arid  differential  measurements  at  TOO  Mev 
correspond  to  typical  balloon  capabilities.  The  values  listed  for 
integrated  intensity  above  30  Mev  refer  to  typical  satellite  measurements. 
Measurements  of  integrated  intensity  above  10  Mev  and  differentia! 
intensity  at  about  30  Mev  represent  the  current  best  satellite  achieve¬ 
ments.  These  are  somewhat  better  than  the  achievements  of  current  radio 
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GENERAL  DETAILS  OF  SOLAR-COSMIC-RAY  EVENTS 

The  number  of  solar-cosmic-ray  events  observed  depends  closely  on 
the  techniques  employed  to  observe  them.  In  Table  2  we  present  a  yearly 
record  of  observed  events  extending  back  to  1942.  Prior  to  1952,  only 
fmj  observations  are  available.  Different  observers  using  this  technique 
agree  quite  closely  on  the  number  of  events  reported.  From  1952  on, 
forward-scatter  observations  are  available.  From  1956,  riometer  obser¬ 
vations  have  been  made.  Both  of  these  meihods  have  yielded  a  yearly 
frequency  of  events  which  is  comparable  to,  but  slightly  less  than,  that 

found  by  f  .  techniques.  Although  the  sensitivities  of  all  these  methods 
min 

are  approximately  the  same,  on  the  basis  of  techniques  alone,  we  may 
falsely  identify  as  solar-cosrnic-ray  events  a  few  instances  of  increased 
ionization  due  to  other  causes. 

Satellite  and  space  probe  observations  have  been  available  since 
1960,  Consequently,  the  number  of  reported  events  has  increased.  The 
total  number  of  events  refers  to  all  separate  events  identified  by  any 
technique.  About  70  events  have  been  observed  at  a  given  sensitivity  of 
detection,  comparable  to  the  sensitivity  of  the  technique,  in  the  solar 
cycle  (cycle  19)  just  being  completed.  In  the  previous  cycle,  approximately 
30  events  were  observed.  On  the  basis  of  measurements  since  1960,  we 
estimate  approximately  100  events  would  have  been  observed  during 
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TABLE  2,  Yearly  Number  of  Solar-Cosmic-Ray  Events  as  Observed 
by  Differerii  Tecl'u liqucs 


r 

Forward 

Year 

min 

r  . 
m  in 

Riometer 

Scatter 

Satellites 

Total 

(Ref. 

17) 

(Ref.  15) 

(Ref.  18) 

(Ref.  16) 

i  (Ref.  19) 

1942 

1943 

0 

0 

1944 

0 

0 

1945 

0 

0 

1946 

6 

6 

1947 

4 

4 

1948 

4 

4 

1949 

•  7 

7 

7 

1950 

3 

5 

5 

1951 

1 

3 

3 

1952 

0 

0 

0 

0 

1953 

0 

1 

1 

1 

1954 

0 

0 

0 

0 

1955 

1 

1 

0 

1 

1956 

5 

6 

2 

3 

6 

1957 

14 

14 

13 

11 

16 

1958 

12 

15 

10 

10 

15 

1959 

7 

8 

5 

5 

1U 

1960 

14 

11 

14 

20 

14*  20** 

1961 

5 

7 

8 

7*  8** 

1962 

1 

3 

1*  3** 

Total 

for  solar 

cycle 

18  (February  1944  -  i 

April  1954)  = 

=  30  events 

Tote  1 

for  solar 

cycle 

19  (April 

1954  - 

H 

'  70  events 
J00  events*** 

*  Not  including  satellite  dofa. 

**  Including  satellite  data. 

***  Estimated  from  satellite  data  1960-1962. 
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cycle  19  to  date,  had  sensitivities  been  the  same  as  in  the  satellite 
space-probe  era.  Although  cycle  19  will  not  end  until  1964,  the  total 
number  of  events  reported  here  are  expected  to  change  very  little,  as  we 
ore  now  In  the  time  of  solar  minimum  when  a  solar  event  is  a  rare  occurrence. 

Of  the  events  reported  during  cyclo  19,  the  52  largest  are  des¬ 
cribed  in  greater  detail  for  the  period  1956-1962  In  Table  3.  It  Is  almost 

certain  that  all  events  with  an  integrated  intensity  greater  than  10^ 

2 

particles/cm  above  30  Mev  are  included  in  the  table  and  events  to  as 

5  2 

low  as  approximately  2x10  particles/cm  above  this  energy  are  included. 
This  is  close  to  the  lower  size  limit,  which  is  presently  determined  by  the 
integral  satellite  techniques.  (This  lower  limit  of  intensity  is  about  equal 
to  the  intensity  of  galactic  cosmic  rays  integrated  over  one  day.) 

Most  of  the  headings  in  the  table  are  self-explanatory,  but  we 
can  make  the  following  remarks  to  elaborate. 

The  integrated  RF  emission  of  a  solar  flare  listed  In  Table  3  is 
obtained  from  intensity-time  profiles.  An  intensity-time  profile  is  o  plot 
of  power  density  per  unit  bandwidth  ct  a  single  frequency,  usually  10,  000 
Mc/s.  The  peak  emission  listed  in  Table  3  is  the  maximum  recorded  on  the 
same  Intensity-time  profile.  The  frequency  of  peak  emission  in  the  microwave 
region  is  not  necessarily  the  same  frequency  at  which  intensity  as  a  function 
of  time  was  recorded. 
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TABLE  3.  Solar  Flare  Data  (1956) 

(Numbers  In  parentheses  refer  to  footnotes  appearing  at  end  of  table) 


Da  te 

February  23 

March  1 1 

August  31 

November  13 

Class 

3+ 

3 

2 

Position^ 

22°N74°W 

n.  d. 

16°N16°E 

16°N10°W 

T  (2) 

opt  max 

0340 

n.d. 

1241 

1301 

Max  RF^ 

20000 
at  0341 

/( RF)(4) 

r  C5) 

'rise 

>  30  Mev 

6 

>  100  Mev 

3 

r  (6) 

'decay 

36(e)°3) 

>  30  Mev 

30 

Uncertain 

Uncertain 

>  1 00  Mev 

16 

Peak  lv7) 

>  10  Mev 

10000 

>  30  Mev 

8000 

150 

>  1 00  Mev 

5000 

60 

Peak  P  ® 
o 

300 

200 

/(l)(9) 

>  10  Mev 

1. 8  x  10 

>  30  Mev 

l.Ox  109 

2.5  x  107 

>  100  Mev 

3.5  x  108 

6x10^ 

<;,  r  ™ 

260 

160 

p/o(,,) 
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TABLE  3.  Sclar  Flare  Data  (1957) 


Date 


Jan  20 


Apr  3 


Ju!  3 


Class 

Position^ 

(2) 

opt  max 

(3) 


T 

Max  RF 


Im{i) 
r.  <5> 

rise 

>  30  Mev 

>1 00  Mev 

r  (6) 

‘ decay 

>  30  Mev 
>100  Mev 

Peak  1^ 

>  10  Mev 

>  30  Mev 
>100  Mev 

Peak  P  ® 

/0)(9) 

>  10  Mev 

>  30  Mev 
>100  Mev 

<K00) 

p/«(u) 


3+ 

25°S30°W 

1120 


24(e) 


2000(e) 

100(e) 

70 


2  x  108 
7  x  106 


3 

154S60°W 

0835 


Uncert 


3+ 

Uncert  14°N40°W 
0740 

2400 


Uncert  30(e) 


200 

50(e) 

2  x  107 


Aug  9 

Uncert 


18(e) 

20(e) 

50(e) 

1. 5  x  106 
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TABLE  3.  Solor  Flare  Data  (1957  Coat'd) 


Date 

Aug  29 

Sep  21 

Oct  20 

Nov  4 

Class 

3 

3+ 

Position^ 

Uncert 

10°N06°W 

26°S45“W 

n.d. 

T  ,  « 

opt  max 

1338 

1642 

Max  RF® 

1100 

/(RF)(4) 

y  (5) 

Mse 

>  30  Mev 

1 

>100  Mev 

r  (6) 
ldecoy 

>  30  Mev 

60(e) 

24(e) 

30(e) 

18(e) 

>100  Mev 

Peak  1^ 

>  10  Mev 

>  30  Mev 

600 

100(e) 

500(e) 

50(e) 

>100  Mev 

30 

140(e) 

Peak  P 

o 

70(e) 

50(e) 

90(e) 

50(e) 

/<o(9) 

>  10  Mev 

>  30  Mev 

1.2  x  108 

1.5  x  106 

5  x  107 

9  x  106 

>100  Mev 

3  x  106 

1  x  107 

‘I*™ 

/  OD 
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TABLE  3.  Solor  Flare  Uara  (1 956} 


Date 

Feb  9 

Mar  23 

Apr  10 

July  7 

Class 

2+ 

3+ 

3+ 

Position^ 

12°S14°W 

14°S78°E 

Uncert 

25®N08°W 

T  (2) 
opt  max 

Max  RF®^ 
/(RF)(4) 

2142 

1005 

10000  at  1004 

0115 

2000  at  01 12 

120 

r  (5) 

1  rise 

>  30  Mev 
>100  Mev 


^decay 


>  30  Mev 

18(e) 

40 

32 

>100  Mev 

20 

16 

Peak  l(?) 

>  10  Mev 

8000  (420) 

1000 

>  30  Mev 

1000(e) 

1200(60) 

50 

1500 

>100  Mev 

100  (1.2) 

100 

Peak  P  (8) 
o 

50(e) 

75 

50(e) 

70 

/( l)(9) 

2  x  109  (9.  5  x  107) 

>  10  Mev 

1.8  x  109 

>  30  Mev 

2.5x  10 8  (7  x  104) 

5x10^ 

2.5  x  108 

>100  Mev 

1  x  107  (8  x  104) , 

9  x  106 

^0 

o 

o 

55 

60 

,  (11) 

p/a 

1.5  ±0.5 

43 


D2-90469 


TABLE  3.  Solar  Flare 

Data  (1958  Cont'd) 

Dote 

Aug  16 

Aug  22 

Aug  26 

Sep  22 

Clots 

3+ 

3 

3 

Position^ 

14*S50#W 

18*N10*W 

20°N54°W 

Uncert 

T  » 

0440 

1448 

0027 

opt  max 

«« i 

Max  RFW' 

6200  at  0442 

6700  at  0025 

/«  F)W 

600 

700 

r.  (i) 

1  rise 

>  30  Mev 

>100  Mev 

10 

10 

3 

r  (6) 

‘decay 

>  30  Mev 

>100  Mev 

18 

20 

8 

12 

30(e) 

Peak  1^ 

>  10  Mev 

3000 

6000 

15000 

600 

>  30  Mev 

400 

600 

1500 

50 

>100  Mev 

20 

20 

50 

1.5 

Peak  P  (3) 

60 

50 

50 

40 

/( D{9) 

A  A  A  — » 

>  10  Mev 

4x  10B 

8  x  10 

1. 5  x  10 

9  x  10 

>  30  Mev 

4  x  107 

7  x  107 

• 

X 

o 

■  CO 

6x10^ 

>100  Mev 

<R(,0) 

1.6  x  106 

1.8  x  106 

2.  Ox  10A 

1  x  105 

40 

/  01) 
p/a 
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TABLE  3.  Solar  Flare  Data  (1959) 


Da  te 

May  10 

June  13 

July  10 

Class 

3+ 

1+ 

3+ 

Position^ 

23°N47°E 

17°N57®E 

22°N70°E 

T  (2) 

2118 

0400 

0222 

opt  max 

Max  RF® 

10000 

20000  at  0224 

/(RF)(4) 

1500 

r.  (5) 

‘rise 

>  30  Mev 

20 

'30 

>1 00  Mev 

10 

24(e) 

18 

^  & 

‘decay 

>  30  Mev 

22 

40 

>1 00  Mev 

14 

20 

Peak  l(7) 

>10  Mev 

30000  (5000) 

15000  (800) 

>  30  Mev 

6000  (500) 

9(e) 

4000  (160) 

>1 00  Mev 

1000  (5) 

1200  (5) 

Peak  P  (8) 
o 

70 

50(e) 

90 

/(l)(9) 

>  10  Mev 

5.5  x  109  (7. 5  x  108) 

4.5  x  109  (1.6  x  108) 

>  30  Mev 

9.6  x  108  (4. 2  x  107) 

8. 5  x  107 

1.0  x  109  (2.  4  x  107) 

>1 00  Mev 

8. 5  x  107  (3, 5  x  105) 

1.4  x  108  (5  x  105) 

70 

75 

p/°(1» 

1.5  ±0.5 

3  ±  1 
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TABLE  3.  Solar  Flare  Data  ('i960) 


Date 

Class 

Position  ^ 

T  (2) 

opt  max 

Max  RF(3) 

/( RF)(4) 

r.  (5) 

‘ti  se 

>  30  Mev 
>100  Mev 

(6) 

^decay 

>  30  Mev 
>1 00  Mev 

Peok  l'71 

>  10  Mev 

>  30  Mev 
>1 00  Mev 

Peak  P  ® 


;o) 


(9) 


Jan  1 1 
3 

23°N03°E 

2100 


«/)P 


30 


2.5 


p/a 


>  10  Mev 

>  30  Mev 

>100  Mev 

00) 

o 

01) 


4  x  10“ 


100 


Apr  1 

Apr  5 

Apr  23 

3 

?<• 

3 

12aN10°W 

10°N61  “W 

05°S34°E 

0859 

0245 

0130 

8000 

1 

3 

,  3  „ 

1 

2 

6 

12 

8 

4 

5 

400 

250 

200 

20 

120 

45 

25 

150 

240 

1,  5  x  107 

1.3  x  107 

i 

5.  Ox  106 

1.1  x  106 

5. 0  x  10 

5 

8.  5  x  10 

7  x  105 

130 

150 
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TABLE  3.  Soiar  Flare 

Da  ia  '  t .  3  Con  I 

1  J\ 

Date 

Apr  29 

May  4 

May 

May  13 

Class 

3+ 

3 

3+ 

Position^  ^ 

Uncert 

14°N90°W 

10oS08°E 

30°N64°W 

T  (2) 
opt  max 

1020 

1450 

0536 

Max  RF® 

200 

8800 

/(RF)(4) 

10 

1200 

r.  (5) 

*  rise 

>  30  Mev 

2-3 

3 

'  3  * 

>100  Mev 

1 

r  (6) 

ldecay 

>  30  Mev 

14 

24 

15 

>100  Mev 

6 

p«k  i<7> 

>  10  Mev 

150 

180 

>  30  Mev 

100 

100 

40 

60 

>100  Mev 

40 

10 

Poak  fj8> 

;o)<9!° 

60 

300 

7 

90 

-y 

>  10  Mev 

1.2  x  10 

1. 5  x  10 

>  30  Mev 

7  x  106 

6  x  106 

4  x  106 

4  x  106 

>100  Mev 

1.2  x  106 

4.  5  x  105 

^5 

o 

o 

200 

85 
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TABLE  3.  Solar  Flare 

Data  (1960  Cont'd) 

Date 

Jun  1 

Aug  12 

Sep  3 

Sep  26 

Class 

3+ 

3+ 

3 

2 

Position^ 

28°N46°E 

22°N27°E 

18°N88°E 

2 1  #S64°W 

T  (2) 

opt  max 

0844 

1940 

0110 

0550 

Mo*  RF® 

12000  at  0108 

h  fim 
r.  (S 

'  rise 

400 

>  30  Mev 

2 

10 

12 

n.d. 

>100  Mev 

6 

r  w 

‘decay 

>  30  Mev 

20 

50 

32 

>100  Mev 

26 

>  10  Mev 

450  (6) 

300 

>  30  Mev 

5 

3 

200  (3) 

40 

>100  Mev 

60  (0.4) 

4 

Peak  P  ® 

0 

180 

65 

/(l)(9> 

>  10  Mev 

9  x  107  (1  x  106) 

2  x  107 

>  30  Mev 

4  x  105 

6x  105 

3. 5  x  107  (3. 6  x  105) 

2  x  106 

>100  Mev 

7  x  106  (4  x  104) 

1. 2  x  105 

<Kao) 

160 

60 

20  ±7 

49 


D2-9Q469 


TABLE  3,  Solar  Flare  Data  (1960  Conf'd) 


Date 

Mev  12 

Nov  15 

Nov  20 

Clou 

3+ 

3+ 

3 

o  W 

Position 

27°N02PW 

30°N32°W 

28°N1 13“W 

T  « 

1329 

0221 

2020 

opt  max 

Max  RF(3) 

10000  at  1329 

14000  at  0227 

1000 

Id rf)(4) 

3400 

2000 

100 

(5) 

Trise 

>  30  Mev 

10 

8 

2 

>100  Mev 

8 

4 

1 

(6) 

Tdecay 

>  30  Mev 

18 

16 

10 

>100  Mev 

14 

10 

5 

Peak  1^ 

22C0 

>  ]0  Mev 

32000  (4000) 

22000  (4200) 

>  30  Mev 

12000  (1500) 

8000  (1500) 

1000 

>1 0C  Mev 

2500  (180) 

2400  (160) 

400 

Peak  P  (8) 
o 

145 

135 

180 

/( ;o(9) 

>  10  Mev 

4  x  109  (4  x  108) 

2.5  x  109  (3  .  8  x  108) 

7.2  x  108  (9  x  107) 

1.2  x  108  (6.  5  x  106) 

1.4  x  108 
7 

>  30  Mev 

1.3  x  109  (1.2  x  108) 

4. 5  x  10 
,„6 

>100  Mev 

2.5  x  108  (1.1  x  107) 

3  x  10 

135 

125 

140 

p/o(,,) 

2  ±0.5 

1.8  i  0.5 
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TABLE  3.  Solar  Flare  Data  (1961) 


Date 

Jul  11 

* 

Jul  12 

Jul  18 

Class 

3 

3+ 

3+ 

Position  ^ 

06°S32°E 

07°S22°£ 

06°S60°W 

T  (2> 

opt  max 

1700 

1030 

1015 

Max  RF® 

2500 

7500 

5000 

/(RF)(4) 

850 

300 

r-  (5) 

‘■rise 

>  30  Mev 

8-10 

8-12 

6-10 

>100  Mev 

4 

24 

2-3 

r  (6) 

‘decay 

>  30  Mev 

22-26 

16-20 

24 

>100  Mev 

18 

12 

12 

Peak 

>  10  Mev 

120 

4000 

7000  (280) 

>  30  Mev 

30 

400 

2500  (100) 

>100  Mev 

3 

15 

60  (11) 

Peak  P  ® 
o 

70 

50 

135 

/(■)(9) 

5  x  108 

>  10  Mev 

1.7  x  107 

1  x  109  (3  x  107) 

>  30  Mev 

3  x  106 

4  x  107 

3  x  108  (5.5  x  106) 

>100  Mev 

2.4  x  105 

1  x  106 

4  x  107  (4  x  105) 

<K(,0) 

60 

40 

110 

6  ±  1 
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TABLE  3.  Solor  Flare  Data  (1961  Cont'd) 


Date 

Jul  20 

Sep  28 

Nov  10 

Class 

3+ 

3 

2+ 

Position^ 

07°S90°W 

U°N30UE 

08°N90°W 

T  (2) 

1600 

2223 

1444 

opt  max 

Max  RF(3) 

/( RF)(4) 

2500 

700 

24 

(5) 

■rise 

>  30  Mev 

4-6 

1.5 

>100  Mev 

1.5 

1 

r  (6) 

•  decay 

>  30  Mev 

6-8 

12 

>1 00  Mev 

4 

8 

Peak  1^ 

>  10  Mev 

400 

500 

>  30  Mev 

150 

150 

>100  Mev 

40 

36 

Peak  P  ® 
o 

150 

210 

/0)(9> 

>  10  Mev 

1.5  x  107 

5  x  107 

>  30  Mev 

5  x  106 

6  x  106 

>100  Mev 

9  x  105 

1. 1  x  106 

</>0,0> 

130 

140 

52 


D2- 90469 


Date 


Class 

Position^ 

T  (2) 

opt  max 

Max  RF^ 
/(RF) (4) 

r.  (5) 

rise 

>  30  Mev 
>100  Mev 


?decoy<6> 

>  30  Mev 
>100  Mev 

Peak 


>  10  Mev 

>  30  Mev 
>1 00  Mev 

Peak  P 

/■(l)(9) 

>  10  Mev 

>  30  Mev 
>100  Mev 


TABLE  3.  Solar  Flare  Data  (1962) 


Feb  4 
3 

10#N80°W  (Uncert) 

0222 


Oct  23 

2 

02#N71°W 

1708 

160 


2 


Uncert 


8 


16 

4.5 


0.5 


80 

6  x  105 
1.2  x  105 
1  x  104 

70 
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TABLE  3.  Solar  Flare  Data  (Cont:d) 


Footnotes 


I 


(1) 

(2) 

a) 


(4) 

(5) 

(6) 

(7) 

(8) 
(9) 

(10) 


Heliographic  position  in  degrees. 
Time  of  optical  maximum  in  UT. 


Maximum  radio  emission  at  a  frequency  in  the  range  from  3000  to 
10000  Mc/s.  The  emission  is  given  in  units  of  10”**  watts/m  -cps 
at  a  time,  if  known.  In  UT. 


Time-integrated  radio  emission  at  a  frequency  jjn  the  range  from 
3000  to  10000  Mc/s  in  units  of  10  joules/m  -cps. 

Onset  plus  rise  time  in  hours  for  two  energy  groups. 

Decay  time  in  hours  for  two  energy  groups. 

Peak  omnidirectional  proton  intensity  in  units  of  particles/cm  -sec  above 
various  energies.  Values  in  parentheses  are  alpha  particle  intensities  in 
units  of  particles/cm*-sec  above  various  energies  per  nucleon. 


Characteristic  peak  P  in  Mv  measured  at  time  of  peak  intensity  of 
10-Mev  particles. 


Time-integrated  omnidirectional  proton  intensity  in  units  of  particles/crr/ 
above  various  energies.  Values  in  parentheses  are  alpha  particle 
intensities  in  units  of  particlos/cm*  above  various  energies  per  nucleon. 


Characteristic  integrated  P  in  Mv  estimated  from  all  available  measure¬ 
ments,  direct  and  indirect. 


(11)  Proton/alpha  ratio  given  for  some  rigidity. 

(12)  No  data. 

(13)  (e)  indicates  estimated  value. 


i 
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1  he  integrated  intensity  above  a  given  energy  is  calculated  using 


fcL  -  ^ ^ ...  j  .1 _ _  .  ,  .. „ _ i r * _ 

me  iulic1  i)i  i c  i  tie  unu  uecu y  lima  vj^cui umy  iy 


J(>  E)  -  /  I  (>  E)  exp(-t/t_)  dt+l  I  (>  E)exp(-t/t  )  dt 
I  max  K  I  max  u 

-Ce  ->o 

■=V'd>U(>e> 

Here  t  =  0  is  the  time  of  maximum  intensity  of  the  appropriate  energy 
particles.  The  values  of  I  are  also  listed  in  the  table.  The  intensity- 
time  characteristics  for  protons  and  alpha  particles  are  assumed  to  be  the 
same.  When  an  event  does  not  show  such  a  characteristically  simple 
intensity-time  behavior,  a  step-by-step  integration  of  the  intensity-time 
profile  has  been  made  to  obtain  the  integrated  intensity  (e.g.,  November 
12,  1960). 

The  integrated  intensities  also  appear  to  exhibit  exponential 

rigidity  spectra.  The  characteristic  Pq  of  this  spectrum  is  close  to,  but 

slightly  smaller  thou,  that  of  the  spectrum  at  peak  intensity  and  is  called 

the  characteristic  integrated  P  . 

o 

When  measurements  are  available  from  two  or  more  techniques,  it 
may  be  possible  to  estimate  roughly  such  quantities  as  decay  and  rise  times, 
peak  intensities,  and  characteristic  rigidities.  These  estimates  are  followed 
by  an  (e)  in  Table  3  and  have  a  greater  uncertainty  than  values  of  these 
quantities  in  events  for  which  more  data  are  available.  In  some  cases,  it 
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is  possioie  to  oscerrain  oniy  rnut  on  cvenl  has  occurred.  Ms  estimates  of 
ilia  can  bo  made  except  tc  say  that  these  event*  ore  generally  smaller  ones. 
We  simply  report  these  as  uncertain. 

A  number  of  features  of  the  solar-cosmic-ray  events  can  be 
deduced  from  Table  3.  Some  of  the  most  interesting  properties  will  be 
described. 


ENORMOUS  VARIATIONS 


An  enormous  variation  in  integrated-particle  intensity  is  observed. 

For  particles  above  30  Mev,  the  integrated  intensity  may  range  from  about 

5  2  9  2 

10  particles/cm  to  about  10  particles/cm.  The  number  of  events  of 

a  given  size,  i.e. ,  the  differential  size  spectrum,  does  not  appear  to 

increase  rapidly  with  decreasing  size  but  seemi  to  peak  at  about  10^ 

2 

particles/cm  and  flatten  off  or  even  decrease  for  sizes  smaller  than  this. 

This  result  may  be  due  in  part  to  a  lock  of  sensitivity  in  detecting  the 

smaller  events  and  in  determining  their  size.  However,  it  may  be  u  real 

effect  rtnee  we  believe  our  probability  of  not  detecting  events  with  an 

6  2 

integrated  Intensity  greater  than  10  particles/cm  at  energies  above  30 

Mev  is  negligible  small.  Furthermore,  an  increase  in  sensitivity,  which 

5  2 

allows  detection  of  events  with  about  10  particles/cm  by  use  of  satellite 
measurements,  causes  about  a  30  percent  increase  in  the  number  of  events 
observed.  We  base  the  30  percent  figure  on  1960-1962  data.  This 
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proportion  of  events,  added  to  the  differential  spectrum  shown  in  Fig.  7, 
mokes  the  distribution  fiui,  bur  is  nor  suggestive  of  any  rapid  increase  in 
number  versus  size  at  smaller  sizes.  This  size  distribution,  shown  also 
in  integral  form  in  Fig.  8,  has  important  implications  fora  statistical 
calculation  of  the  integrated  intensities  of  particles  received  over 
relatively  long  periods.  In  effect,  such  a  distribution  implies  that  the 
probability  of  receiving  a  dangerous  dose  depends  on  the  probability  of 
encounter  with  a  large  event  and  that  the  addition  of  small  events  Is 
not  a  serious  problem. 

ASSOCIATION  WITH  ACTIVE  CENTERS 

We  note  that,  although  52  separate  events  are  catalogued,  they 

occurred  in  only  34  active  centers.  During  the  three  years  from  1959 

through  1961,  this  effect  was  even  more  pronounced:  27  sepa rave  events 

are  associated  with  15  active  centers.  Usually,  the  active  centers 

associated  with  the  cosmic-ray-producing  flares  were  already  well 

developed  when  they  appeared  on  the  east  limb  of  the  sun.  Only  once 

after  a  major  event  has  the  same  active  center  continued  to  produce 

cosmic-ray  events  on  its  subsequent  appearance  on  the  visible  portion 

of  the  sun.  We  are  thus  motivated  to  group  the  cosmic-ray  events 

according  to  active  centers.  During  solar  cycle  19,  eight  centers  have 

9  2 

each  produced,  within  o  period  of  a  week,  in  excess  of  10  porticles/cm 
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above  lOMev  ot  the  earth.  The  intensities  produced  by  these  centers 
exceed  the  yearly  intensities  of  galactic  cosmic  rays  by  a  factor  of  ten. 

The  output  of  cosmic  rays  above  10  Mev  from  these  eight  active  centers 
is  about  90  percent  of  all  the  solar  cosmic  rays  above  this  energy  produced 
during  this  cycle.  These  eight  major  events  are  listed  year  by  year  in 
Table  4.  Although  none  of  these  major  events  occurred  in  1957  at  the 
peak  of  the  sunspot  cycle,  the  correlation  of  the  frequency  of  these 
events  with  sunspot  number  is  much  better  than  the  similar  correlation  of 
the  frequency  of  events  observed  at  sea  level.  Sec-level  events  appear 
to  occur  more  frequently  on  the  increasing  and  decreasing  parts  of  the 
cycle  (Ref,  2). 

EAST-WEST  ASYMMETRY 

Of  events  for  which  we  are  able  to  make  a  size  estimate  and  also 
to  know  the  location  of  the  producing  flare  on  the  sun,  29  have  come  from 
flares  on  the  sun's  western  hemisphere  and  8  from  the  eastern.  Of  these 
same  events,  24  have  come  from  flares  in  the  northern  hemisphere  and  12 
from  the  southern.  The  north-south  asymmetry  aepenas  partly  on  a  few 
active  centers  and  may  not  be  real.  The  east-west  asymmetry,  which  has 
been  pointed  out  previously  (Refs.  3  and  6),  is  of  great  significance, 
however.  This  effect,  which  is  undoubtedly  a  propagation  effect  through 
interplanetary  space,  holds  for  all  size  groups,  not  just  the  largest  events. 
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TABLE  4. 

The  Eight  Largest  Solar-Cos 

mic-Ray  Outbursts  During  Cycle  19 

Year 

Number  of  Events 

L>afes 

1954 

0 

1955 

0 

1956 

1 

February  23 

1957 

0 

1V58 

3 

March  26,  July  7,  August  16-26 

1959 

2 

May  10,  July  10-16 

1960 

1 

November  12-15 

1961 

1 

July  12-18 

1962 

0 

1963 

0 
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This  peculiar  result,  which  may  follow  from  the  rather  flot  size  distribution 
of  events,  means  that  flares  occurring  in  active  regions  on  the  eastern 
hemisphere  of  the  sun  are  about  one-third  as  likely  to  produce  a  cosmic 
ray  event  (of  any  size)  observed  at  the  earth.  In  addition,  the  average 
onset-delay  and  rise  times  for  cosmic  rays  from  flares  on  the  east  hemisphere 
of  the  sun  are  about  three  times  as  long  as  from  the  western  hemisphere. 

This  situation  suggests  that,  on  the  average,  a  greater  fraction  of 
the  emitted  particles  is  received  at  the  earth  when  the  flare  is  on  the 
western  hemisphere  of  the  sun.  - 

ALPHA  PARTICLES  AND  HEAVIER  NUCLEI 

In  many  events,  alpha  particles  and  heavier  nuclei,  principally 
C,  N,  O,  and  Ne,  hove  been  observed  in  considerable  numbers  (Refs.  7 
and  8).  These  particles  appear  to  have  the  same  rigidity  spectrum  as  the 
protons  and  have  essentially  the  same  intensity-time  characteristics  as  the 
protons  if  the  comparison  is  on  a  rigidity  basis.  These  facts  suggest  that 
the  interplanetary  control  mechanism  sorts  the  particles  according  to 
rigidity.  The  ratio  of  protons  io  alphas  above  a  fixed  rigidity  varies 
widely,  whereas  the  alpha  to  heavier  nuceii  ratio,  based  on  only  three 
events,  appears  to  be  a  constant,  approximately  sixty.  Most  commonly, 
the  value  for  the  proton-alpha  ratio  lies  between  one  and  two,  although 

events  with  a  ratio  around  20  have  been  observed.  Even  at  a  ratio  of 
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:  s.!  !.<  ui'lev  I  :.ini_L\  ui  a  fixed  rig  laity, 

•  i  •.  ••r.i.o  "«:■!'/  less  i hot  of  protons.  This  detection  problem 
particiJoriy  hampers  balloon  observations,  and,  for  this  reason,  the  alphas 
have  been  observed  only  during  the  largest  events.  It  is  quite  possible 
that  the  proton-alpha  ratio  increases  with  decreasing  event  size.  Then 
the  preponderance  of  ratios  in  the  range  from  one  to  two  would  be  mis¬ 
leading.  In  fact,  there  is  some  evidence  for  this  increase.  In  at  least 
two  of  the  smaller  events  (July  12,  1961  and  September  28,  1961),  an 
absence  of  observed  alphas  implies,  within  the  limits  of  detection,  proton- 

alpha  ratios  greater  than  five.  Integrated  fluxes  of  alphas  with  energy 

8  2 

greater  than  30  Mev/nucleon  were  observed  to  be  about  2. 5  x  10  /cm 

during  the  July  1959  events.  During  the  November  1960  events,  this 

8  2 

intensity  reached  2.  2  x  10  /cm  .  These  values  amount  to  about  10  per¬ 
cent  of  the  flux  of  protons  with  the  same  energy  range  in  these  events 
and  to  about  three  times  the  total  galactic-particle  flux  integrated  over 


one  year.  In  the  November  1960  event,  we  assume  that  the  particles 


heavier  than  alphas  h'-r‘  the  same  intensity-time  profile  as  the  alphas. 


Then  the  total  integrated  intensity  of  such  nuclei  with  energies  greater 

6  2 

than  30  Mev/nucleon  was  about  3x10  /cm  . 
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SUMMARY  OF  OBSL  RVrD  iKTtNSI  TIES  OF  SOLAR  COSMIC  RAYS 


By  way  of  a  general  summary  of  the  intensities  of  solar  cosmic  rays, 

we  show  in  Table  5  the  yearly  integrated  intensities  of  solar  protons  above 

10,  30,  and  100  Mev,  respectively.  A  similar  table  could  be  made  for 

alphas,  but  since  our  information  relates  to  only  a  few  events,  we  have 

not  attempted  to  make  this  compilation.  It  is  seen  from  the  table  that  the 

peak  intensity  of  protons  was  observed  in  1959  for  all  energy  intervals. 

This  maximum  is  due  principally  to  the  July  events  and  again  points  up 

the  importance  of  the  single  large  event.  The  intensity,  integrated  over 

the  entire  solar  cyclo  19,  is  very  nearly  equal  to  the  sum  over  the  six  years 

9  2 

shown  in  the  table  and  amounts  to  9  x  10  particles/cm  above  30  Mev  in 

11  .  2 

free  space.  We  estimate  that,  over  the  entire  cycle,  2  x  10  particles/cm 

above  3  Mev  would  be  observed  near  the  earth.  Also  shown  in  this  table 

are  the  P  's  characteristic  of  the  spectra,  based  ori  the  yearly  integrated 
o 

intensities.  Notice  that  the  P  's  are  smaller  (50  to  70  Mv)  during  the 

o 

years  near  sunspot  maximum.  Accordingly,  no  sea-level  events  are 

observed  during  these  years.  However,  the  total  yearly  flux  of  oil  solar 

cosmic  rays  is  highest  during  these  years.  The  flatter  spectra  observed 

during  the  ascending  and  descending  portions  of  the  solar  cycle  are  in 

keeping  with  observation  of  the  greatest  frequency  of  sea-level  events 

during  these  periods,  although  the  total  yearly  intensity  is  then  somewhat 

64 


TABLE  6.  Observed  Yearly  Inteijio  led  Intensities  at  '»uiui  IVituiis 


Integrated  Iniensi 

11 _ 

Year 

Number  of  Events 

E  >  10  Mcv 

E  >  30  Mcv 

E  >  1  00  Mev 

P 

O 

(protons/cm^) 

(Mv) 

1956 

4 

2.0  x  109 

1.0  x  109 

3.5  x  108 

250 

1957 

9 

— - 

4.0  x  108 

2.  Ox  107 

60 

1958 

8 

7.  Ox  109 

7.  8  x  108 

2.  4  x  107 

50 

1959 

6 

2. 2  x  1010 

4. 2  x  109 

4. 6  x  1 08 

70 

1960 

15 

6.8  x  109 

9 

2.2  x  10 

3.8  x  108 

130 

1961 

6 

1. 6  x  109 

3. 5  x  108 

4. 2  x  107 

90 

1962 

2 

1963 

_J 

52 

Solar-cycle  total 

3. 9  x  1010 

9.  Ox  109 

1.3  x  109 

100  (ay. ) 
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less  fiiOi i  iii  t/iu  peak  period. 

foi  environments  noOr  the  surface  of  the  earth,  these  oulj  are 

presented  in  another  way  in  Fig.  9.  Here,  the  intensity  is  plotted  as  a 

function  of  geomagnetic  latitude  and  the  parameter  L,  which  are  related 

2 

by  cos  X  =  1/L. 

In  Fig.  9  the  curve  labelled  (a)  is  the  integrated  intensity  of 
particles  of  all  energies  during  the  "worst"  year  in  the  recent  sunspot 
cycle  (flux  in  1956  at  high  energies;  1969  at  low  energies);  (b)  is  the 
average  yearly  intensity  observed  during  the  recent  sunspot  cycle; 

(1)  is  the  envelope  of  the  greatest  single  events,  consisting  of  the  July 
1959  events  at  low  energies  (<  50  Mev),  the  November  1960  events  at 
intermediate  energies  (50  to  150  Mev),  and  the  February  1956  event  at 
higher  energies  (>  150  Mev);  (2)  is  the  envelope  of  the  second-largest 
single  events,  consisting  of  the  November  1960  events  at  low  energies, 
the  July  1959  events  at  intermediate  energies,  and  the  November  1960 
events  again  at  higher  energies;  and  (3)  is  the  envelope  of  the  third- 
largest  single  events,  consisting  of  the  May  1959  event  at  low  energies, 
February  1956  at  intermediate  energies,  and  the  July  1959  events  at 
higher  energies. 

These  plots  show  clearly  the  importance  of  the  single  large  events 
as  well  as  the  effects  of  the  particle  spectra  in  the  different  events. 
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Finally,  we  show  in  Fig.  9  the  yearly  integrated  intensities  of 

all  galactic  cosmic  rays  to  be  expected  at  both  sunspot  minimum  and 

sunspot  maximum.  We  see  that  above  a  nominal  geomagnetic  latitude  of 

65°  or  outside  a  magnetic  shell  with  L  =  5.  5  (i.  e. ,  in  the  equatorial  plane 

5.5  times  the  radius  of  the  earth),  the  intensities  of  solar  cosmic  rays 

dominate  those  of  galactic  particles  in  yearly  averages.  This  domination 

holds  true  in  all  but  the  quietest  years  near  sunspot  minimum. 

An  important  point  should  now  be  mentioned  in  connection  with 

the  spatial  distribution  of  solar  cosmic  rays  in  the  earth's  magnetosphere. 

Experimental  flux  measurements  are  shown  as  a  function  of  the  conventional 

geomagnetic  latitude  on  the  lower  scale  of  Fig.  9.  From  these  measurements 

and  with  certain  assumptions  of  isotropy,  an  energy  or  rigidity  spectrum  of 

the  cosmic  rays  could  be  constructed  if  the  geomagnetic  field  were  at  all 

times  a  simple  dipole  derivable  entirely  from  sources  within  the  earth. 

For  example,  if  the  geomagnetic  field  is  assumed  to  be  a  dipole  with  a 
25  3 

moment  of  8. 1  x  10  gauss-cm  ,  the  cutoff  (minimum)  rigidity,  P  ,  for 

c 

particles  reaching  the  earth's  surface  at  a  given  latitude  X  is  given  by 

P  =  14.  9  cos^  X  Gv.  v 
c 

Such  calculations  agree  reasonably  well,  at  times  when  the  sun  is  quiet 
(i.  e. ,  no  solar  cosmic  rays),  with  direct  spectsa  I  measurements.  However, 
during  solar  event).,  ring  currents  in  the  magnetosphere  and  magnetic  storms 


duo  to  currents  just  outside  Inc  mcj; »i is.-1  o--; »i i * -i  e  tfuciuuse  rile  cutott 
rigidity  at  a  given  point.  Calculations  of  the  spectrum  using  geomagnetic 
coordinates  would  then  be  incorrect.  On  the  average,  during  the  arrival 
of  solar  cosmic  rays  at  the  earth,  the  cutoff  rigidities  at  latitudes  greater 
than  60°  may  be  only  about  one-half  of  what  is  normally  to  be  expected. 
Thus,  detecting  stations  receive  particles  which,  on  the  basis  of  dipole 
field  calculations,  would  be  expected  only  at  latitudes  about  5°  higher 
than  the  actual  geomagnetic  latitude  of  the  station.  Equivalently,  the 
solar  particles  reach  lower  latitudes  and  smaller  L  values  than  normally 
expected.  This  effect  is  variable  from  event  to  event,  but  an  average 
correction  is  deemed  suitable  for  the  preceding  kind  of  analysis  and  is 
shown  as  a  set  of  effective  latitudes  and  L  values  at  the  top  of  Fig.  9. 

The  net  effect  is  to  move  the  latitude  at  which  the  solar  cosmic  rays 
dominate  the  galactic  particles  to  about  60°  as  measured  in  the  actual 
geomagnetic  coordinates  instead  of  65°  as  predicted  theoretically,  and  to 
change  L  ftom  approximately  5.5  to  4,  As  a  consequence,  the  magnetically- 
shielded  region  around  the  earth  is  shrunk. 
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PREDICTION  CAPABlLlllhi  (-Ok  SOLAR-CCSMiC-RAY  EVENTS 

The  ability  to  predict  the  occurrence  of  a  solor-cosmic-ray  event 
is  of  prime  importance  in  planning  space  missions.  Prediction  capabilities 
may  be  resolved  into  three  rather  distinct  categories. 

FEW-HOUR  WARNING 

Solar  cosmic  rays  generally  do  not  reach  peak  intensity  until  a 
number  of  hours  after  the  source  flare.  If  we  can  estimate  reasonably  from 
the  size  of  a  flare  or  from  its  accompanying  electromagnetic  emissions  the 
integrated  cosmic-ray  intensity  to  arrive  at  the  earth,  then  we  could  pro¬ 
vide  a  warning  several  hours  in  advance  of  the  solar-cosmic- ray  event. 

During  this  time,  astronauts  could  make  special  preparations  and  perhaps 
move  an  orbiting  vehicle  to  less  dangerous  regions  of  space. 

If  we  attempt  to  use  flare  size  as  a  criterion,  'hen  we  have  the 
following  data  from  solar  cycle  1?.  During  this  cycle,  a  total  of  about 
120  flares  of  class  3  or  larger  was  observed.  Of  these,  33  were  accompanied 
by  solar-cosmic-ray  outbursts.  In  this  same  period,  of  the  approximately 
1,  100  class  2  flares  which  occurred,  6  were  accompanied  by  cosmic-ray  out¬ 
bursts.  Obviously,  there  is  some  correlation  with  flare  size,  but  it  is  not 
a  very  strong  criterion.  We  need  an  indicator  that  will  give  not  only  a 
high  probability  of  predicting  the  occurrence  of  a  eosmic-roy  event,  but 
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also  some  measure  of  the  size  of  the  cosmic-ray  event.  We  believe  that 
such  a  measure  is  provided  by  the  correlation  between  the  microwave 
(cm  wave)  output  from  the  flare  and  the  cosmic-ray  intensity  observed 
at  the  earth,  in  particular,  we  shall  use  for  our  correlation  study  the 
emission  at  10,  000  Mc/s.  In  Fig.  10,  we  show  a  regression  curve  of  the 
peak  emission  at  10,000  Mc/s  and  the  peak  intensity  of  cosmic  rays  above 
10  Mev  at  the  earth.  We  see  a  strong  correlation  between  these  quantities; 
however,  some  events  lie  unusually  far  off  the  curve.  A  number  of  features 
could  cause  this  scatter.  First  of  all,  this  approach  does  not  consider  the 
obvious  spectral  differences  of  the  events.  Hence,  comparison  of  radio 
emission  at  a  fixed  frequency  and  intensity  of  particles  above  a  given 
energy  may  give  different  relations  for  different  spectra.  A  regression 
curve  similar  to  that  in  Fig.  10,  but  using  the  peak  intensity  of  cosmic  rays 
above  30  Mev,  shows  a  wider  scatter. 

In  addition,  the  intensity  of  cosmic  rays  near  the  earth  depends 
strongly  on  the  propagation  conditions  in  interplanetary  space  and  on  the 
location  of  the  flare  on  the  sun.  An  east-limb  flare  might  be  expected  to 
produce  less  cosmic-ray  intensity  at  the  earth  than  a  west-limb  flare  of 
equal  magnitude,  for  example. 

|p  view  of  these  difficulties,  which  are  inherent  in  any  warning 
process  based  in  this  general  way  on  electromagnetic  emission,  we  have 
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compared  in  Fig.  1 1  the  integrated  emission  at  10,000  Mc/s  with  the 
integrated  intensity  of  cosmic  rays  above  10  Mev  observed  at  the  earth. 

The  correlation  here  is  perhaps  somewhat  better  than  for  the  peak  intensities, 
and,  of  course,  the  integrated  intensity  is  a  more  useful  quantity  for  the 
dose  calculations.  In  fact,  if  we  omit  the  two  events  which  give  one  to 
two  orders  of  magnitude  less  cosmic-ray  Intensity  than  the  prediction 
Indicated  by  their  radio  emission,  we  find  that  we  can  predict  within  a 
factor  of  three  the  integrated  cosmic-ray  intensity,  given  the  integrated 
radio  emission. 

PREDICTIONS  OF  A  FEW  WEEKS  OR  LESS 

The  capability  to  predict  a  few  weeks  or  less  in  advance  would  be 
valuable  for  orbital  or  lunar  missions  of  this  duration.  We  would  use 
visible  conditions  on  the  sun  and  the  possible  occurrence  of  a  cosmic-ray- 
producing  event  as  basic  factors.  As  we  have  seen,  the  presence  and 
development  of  an  active  region  with  its  associated  sunspots  and  complex 
magnetic  field  is  a  basic  port  of  the  process  which  leads  to  a  solar-cosmic- 
ray  event.  We  have  also  seen  that  the  larger  events  constitute  the 
greatest  hazard  for  any  mission,  at  least  for  one  of  extended  duration.  Let 
us  first  examine  the  histories  of  twelve  aciive  centers  which  went  on  to  pro¬ 
duce  the  largest  cosmic-ray  events.  Of  these  twelve  centers,  three 

proauced  major  cosmic-ray  events  on  their  first  passage  across  the  visible 
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solar  hemisphere,  seven  on  their  second  passage,  and  one  each  on  the  third 
and  fifth  recorded  passages.  The  average  time  from  the  birth  of  an  active 
center  to  the  production  of  solar  cosmic  rays  thus  seems  to  be  a  little  less 
than  two  solar  rotations  (54  days). 

Those  events  occurring  on  the  first  passage  can  give,  at  most,  a 
few  days'  warning.  In  all  these  instances,  the  active  center  was  fairly 
large  and  growing  in  size  and  intensity  as  it  appeared  at  the  east  limb,  in 
five  cases  where  events  occurred  on  the  second  passage,  the  associated 
active  centers  were  sufficiently  large  and  growing  to  seem  as  potentially 
dangerous  as  the  centers  which  produced  events  on  their  first  passage. 

The  other  two  centers  which  produced  second-passage  events,  including 
the  November  1960  major-event  group,  were  relatively  inconspicuous.  Of 
the  two  events  occurring  on  the  third  and  fifth  passages  of  the  associated 
active  centers,  the  one  in  July  1959  was  associated  with  a  large  and 
growing  center  but  the  other.  May  1959,  was  associated  with  a  weak 
center  which  persisted  for  four  rotations  and  gave  no  special  indications 
before  it  burst  into  cosmic-ray  activity  on  its  fifth  rotation. 

In  addition,  only  two  active  centers  have  produced  cosmic  rays 
on  more  than  one  rotation. 

All  of  the  eight  largest  events  in  cycle  19  occurred  on  the  second 
or  a  later  passage  of  the  associated  active  center.  This  gives  us  at  least 
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q  possibility  of  recognizing  the  potentially  dangerous  event*  perhaps  a 
month  ahead  of  time.  Unfortunately,  no  reliable  method  has  been  devised 
for  telling  which  of  these  active  regions  is  potentially  dangerous.  We 
must  resort  to  the  simplest  possible  criterion,  a  criterion  based  on  the  size 
of  the  active  center  or  of  the  sunspots  making  up  the  active  center.  If 
we  adopt  such  an  approach  and  set  our  criterion  too  low,  then  it  is 
possible  that  no  safe  periods  at  all  will  be  found  near  sunspot  maximum. 
Application  of  this  approach  to  selected  years  by  Anderson  (Ref.  20) 
shows  that  adopting  a  critical  area  of  1000  x  10  ^  of  that  of  the  sun  reduces 
the  usable  time  to  '**  25  percent  near  sunspot  maximum.  This  criterion  also 
reduces  the  number  of  events  encountered,  but,  unfortunately,  not  by  a 
great  deal  more  than  the  factor  of  four  reduction  in  time.  A  study  of 
selection  criteria  of  this  sort  for  various  area  limits  reveals  that  if  we  are 
wi  lling  to  accept  a  usable  time  between  25  and  50  percent  of  that  avail¬ 
able,  then  we  can  achieve  over  the  sunspot  cycle  a  factor  of  two  to 
three  decrease  in  frequency  of  solar-cosmic-ray  events  by  means  of  the 
size  criterion. 

Another  method  of  decreasing  the  frequency  of  encounter  with  a 
solar-cosmic-ray  event  makes  use  of  the  east-west  asymmetry  in  the 
location  of  the  flare  on  the  sun  and  the  cosmic-ray  events  observed  at 
the  earth.  Figure  12  shows  the  number  of  cosmic-ray  events  observed  at 
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FIGURE  12.  Occurrence  of  solar-cosmic -ray  events  c  the  earth  as  a 
function  of  longitude  of  the  flare  on  the  visible  solar 
hemisphere. 
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the  earth  versus  the  longitude  on  ihe  sun  of  the  producing  flare  foi 
cycle  19.  We  see  that  the  probability  of  seeing  cosmic  rays  from  flares 
or  active  centers  on  the  west  quadrant  is  about  three  times  greater  than 
that  from  flares  on  the  east  quadrant.  Missions  may  thus  be  conducted 
during  the  period  where  an  active  center  is  on  the  east  quadrant  with  a 
decreased  probability  of  encounter  with  a  cosmic-ray  event. 

In  summary,  we  may  say  that  without  a  great  deal  of  understanding 
of  the  processes  causing  cosmic-ray  flares,  we  may  still  decrease  our 
probability  of  encounter  by  a  factor  of  five  for  intervals  of  one-  to  two- 
weeks'  duration  at  the  expense  of  reducing  the  available  time  by  50  percent. 
The  necessary  predictions  can  be  made  perhaps  30  days  in  advance  by 
using  relatively  simple  criteria  such  as  development  time,  location,  and 
sire  of  flare-producing  active  centers. 

Other  systematic  studies  of  the  active  centers  may  improve  this 
prediction  capability  somewhat,  but  our  understanding  of  the  basic  solar 
physics  must  be  improved  to  provide  a  rigorous  basis. 

PREDICTIONS  FROM  MONTHS  TO  YEARS 

A  capability  to  predict  for  time  intervals  in  excess  of  a  few  months 
and  extending  up  to  a  number  of  years  is  needed  to  plan  missions  with  lunar 
bases  and  for  planetary  missions  which  may  take  more  than  a  year.  This 
capability  is  necessarily  based  on  an  extrapolation  of  our  knowledge  of 
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the  variation  in  the  rate  of  occurrence  of  solar-cosmic-ray  events  during 
the  sunspot  cycle. 

Let  us  first  examine  some  relations  between  the  solar-cosmic-ray 
events  and  certain  features  relating  to  the  sunspot  cycle.  There  are  a 
number  of  solar  features  we  could  examine.  These  include  (1)  the  sunspot 
number,  (2)  the  yearly  distribution  of  3+  flares,  and  (3)  the  yearly  distribu¬ 
tion  of  microwave  bursts  above  a  certain  intensity.  Even  though  the  sun¬ 
spot  number  is  not  a  good  indicator  of  individual  events,  we  suggest  that 
sunspot  number  may  be  a  useful  indicator  of  the  yearly  integrated  intensify 
of  solar  cosmic  rays  above,  say,  10  Mev.  The  satisfactory  correlation 
between  these  quantities  is  shown  in  Fig.  13  for  the  years  from  1955 
through  1962.  This  indicates  that  we  can  make  estimates  of  the  expected 
yearly  intensities  of  cosmic  rays  from  projections  of  yearly  sunspot  numbers. 
To  investigate  this  correlation  further  and  to  extend  the  measurements  over 
more  than  one  solar  cycle,  we  show  in  Fig.  14  a  plot  of  the  yearly  number 
of  solar-cosmic-ray  events  as  a  function  of  the  yearly  sunspot  numbers  for 
the  years  from  1942  through  1962.  The  number  of  cosmic-ray  events  is 
based  on  f  .  data  adjusted  to  a  constant  sensitivity  of  detection  over  the 
entire  period.  As  we  have  pointed  out,  this  adjustment  should  amount  to 
at  least  a  30  percent  increase  to  make  the  data  comparable  with  that  pro¬ 
vided  in  the  era  of  improved  detection  capabilities  provided  by  satellites 
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INTEGRATED  INTENSITY  OF  PARTICLES  ABOVE  10  Mev  (particles/cm2) 


FIGURE  13.  Smoothed  sunspo'  number  as  a  function  of  integrated  intensity  of 
solar  cosmic  rays  above  10  Mev. 
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FIGURE  14.  Number  of  solar-cosmic-raj'  events  observed  as  a  function  of 
the  smoothed  sunspot  number. 
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and  space  probes 

As  can  be  seen  from  Fig.  14,  there  is  c  usable  degree  ot  linear 

correlation  between  the  two  parameters.  In  particular,  we  see  that 

during  ten  years  when  the  yearly  sunspot  number  has  been  less  than  forty, 

only  three  solar-cosmic-ray  events  have  been  observed — an  average  of 

0.3  events/year.  Further  from  Fig.  13,  we  see  that  the  total  yearly 

integrated  intensity  of  particles  above  10  Mev  expected  from  these  events 

8  2 

is  less  than  approximately  10  /cm  ,  or  approximately  the  integrated 
intensity  of  galactic  cosmic  rays  in  one  year. 

Observations  over  the  past  three  sunspot  minima  show  that  the 
sunspot  number  remains  below  forty  for  an  average  of  three  years  around 
each  minimum.  During  the  current  period  of  sunspot  minimum,  which 
started  early  in  1962  and  will  probably  extend  to  early  in  1965,  travel 
outside  the  earth's  magnetosphere  should  be  relatively  safe. 

The  character  of  the  next  sunspot  maximum  is  the  subject  of  some 
study  in  the  literature.  According  to  the  work  of  Gleissberg  (Ref.  21), 
the  next  maximum  will  occur  late  in  1967  or  early  in  1968,  and  it  will 
have  a  number  of  approximately  85.  We  recognize  that  the  sunspot  number 
of  cycle  19  was  the  largest  on  record;  in  fact,  the  maximum  sunspot  number 
during  cycle  18  was  one  of  the  largest  on  record.  Thus,  apparently 
1 1-year  cycles  of  high  sunspot  numbers  are  coming  to  an  end  in  accordance 
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with  the  80-ycar  sunspot  cycle.  If  we  can  extrapolate  the  correlations 

observed  in  cycles  18  and  19  to  cycle  20  and  accept  Gleissberg's  estimates 

of  maximum  sunspot  number,  in  1967  and  1968  we  might  expect  five  or 

six  soiar-cosmlc-ray  events  per  year  and  a  yearly  integrated  inter. ity  of 

9  2 

particles  above  10  Mev  of  about  10  /cm  .  Although  the  probability  of 
obtaining  events  with  flatter  spectra  appears  greater  due  to  the  lower 
sunspot  number,  generally  the  hazard  during  the  next  period  of  sunspot 
maximum  appears  to  be  substantially  reduced  from  the  hazard  during  the 
previous  maximum. 
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SOLAR-PARTICLE  EVENT  RADIATION  DOSE  DETERMINATION 

Previous  sections  of  this  report  have  established  that  the  principal 
radiation  threat  in  the  polar  regions  and  the  regions  beyond  the  magneto¬ 
sphere  is  the  solar-particle  event.  In  order  to  determine  the  magnitude  of 
this  threat  for  missions  in  these  regions  of  space,  it  is  necessary  to  calculate 
the  radiation  doses  encountered  within  presently  designed  spacecraft. 

This  section  of  the  report  is  devoted  to  a  description  of  the  method  used  in 
calculating  the  radiation  doses  for  solar-particle  events. 

in  general,  the  method  employed  is  that  developed  by  D.  L.  Dye 
(Ref.  22).  In  this  method,  the  unit  flux  dose  or  dose  per  proton  is  calculated 
by  using  a  computer  program.  With  this  computer  program,  the  unit  flux 
dose  may  be  determined  at  any  point  in  any  material  /or  which  a  value  of 
dE/dx  Is  available.  This  material  may  be  in  free  space  or  shielded  by  as 
many  as  twelve  layers  of  different  materials.  In  addition,  the  tissue  depth 
dose  may  be  determined  by  using  a  subroutine  employing  a  body  phantom. 

The  unit  flux  dose  includes  primary  protons  and  secondary  radiations 
generated  in  the  shielding  materials  and  the  body  tissue.  The  radiation  for 
the  solar  particle  event  is  then  determined  by  multiplying  the  unit  flux  dose 
times  proton  flux  for  the  event. 

The  first  step  in  calculating  the  unit  flux  dose  is  the  specification 

of  particle  spectrum.  From  previous  sections,  it  has  been  shown  that  the 
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radiation  of  the  solar  particle  event  is  predominantly  protons  with 
energies  from  one  Mev  to  a  few  Gev.  The  integral  spectrum  best 

describing  the  flux  for  a  given  event  is: 

J(>P)  *  Jq  exp (-y-), 


(Eq.  1) 


where  P  is  the  particle  rigidity  ar>d  J(>  P)  is  the  number 

of  particles  of  rigidity  greater  than  P.  The  e-folding  rigidity  Pq  varies 
from  event  to  event,  ranging  from  50  to  200  Mv.  In  Table  6,  the  estimated 
integral  flux  at  10,  30,  and  100  Mev  is  listed  for  those  events  for  which 
sufficient  data  are  available.  For  some  event.',  only  an  estimate  of  the 
flux  above  30  Mev  can  be  made.  In  order  to  proceed,  a  value  of  Pq 
for  these  events  must  be  determined.  This  is  accomplished  as  follows. 


Substitution  of 


p=-  y?+ 

q 


(Eq.  2) 


into  Eq.  1,  where  T  =  kinetic  energy,  q  =  charge,  and  m  =  rest  mass  in 
energy  units,  gives  the  spectrum  as  a  function  of  T,  i.  e. , 


J'(>  T)  =  Jo  exp 


+  2m  T  . 


(Eq..  3) 


Substitution  of'T  =  30  Mev  and  T  =  100  Mev  into  Eq.  3  gives,  for  protons: 

J»-|>foO)  “  ^P  F  10°2+  187600  "  3°2  +  56280  * 

(Eq.  4) 
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TABLE  6.  Integral  Flux  at  10,  30,  and  100  Mev  and  Corresponding 
Characteristic  Kigidity  F 

o 


Date 

J(>  10  Mev) 

J(>  30  Mev) 

J(>  100  Mev) 

P 

0 

2/23/56 

1. 8  x  109 

(protons/cm^) 
1.0  x  109 

3.5  x  108 

(Mv) 

195 

3/1 1/56 

- 

7 

- 

8/  3/56 

- 

2.5  x  10 

6  x  10 

144 

11/13/56 

- 

- 

- 

- 

1/20/57 

- 

2  x  108 

7  x  10° 

61 

4/  3/57 

- 

- 

- 

- 

6/22/57 

- 

7 

- 

- 

7/  3/57 

- 

2x1°; 

- 

- 

8/  9/57 

- 

1. 5  x  10° 

1.2  x  10° 

- 

8/29/57 

- 

3  x  10° 

56 

9/21/57 

- 

1.5  x  10° 

1  x  107 

- 

10/20/57 

- 

5xi°; 

127 

11/  4/57 

- 

9  x  10 

- 

- 

2/  9/58 

-  o 

1  X  l0fi 

2. 5  x  10° 

7 

- 

3/23/58 

2  x  10V 

1  X  10 

64 

4/10/58 

o 

5x  10° 

2.5  X  10T 

A 

- 

7/  7/58 

1. 8  x  10* 

4x  10° 
8x10° 

9x10° 

62 

8/16/58 

4  x  lo' 

1.6  x  10° 

64 

8/22/58 

7  x  lo' 

1.1  x  10° 

1.8  x  10° 

56 

8/26/58 

1. 5  x  10* 

2.  Ox  10° 

51 

9/22/58 

9  x  10 

6x!06 

1  x  10 

50 

5/10/59 

5. 5  x  109 

9. 6  x  i  0? 

8. 5  x  107 

84 

6/13/59 

*"  Q 

8.5  x  10' 

1.4  x  108 

1.0  x  10° 

1.3  x  10° 

- 

7/10/59 

4. 5  x  1  Og 

1.0  x  10* 

104 

7/14/59 

7. 5  x  10* 

1.3  x  10* 

9.1  x  10° 

80 

7/16/59 

3.3  x  10 

105 

8/18/59 

- 

1.8  x  10° 

- 

- 

1/1 1/60 

7 

4x,0f 

5 

- 

4/  1/60 

1.5  x  10 

5.  Ox  10° 

8.5  x  10° 

116 

4/  5/60 

7 

1.1  x  10° 

5 

- 

4/28/60 

1.3  x  10 

5.0x10° 

7  x  10 

104 

4/29/60 

-  • 

7  x  10° 

- 

- 
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TABLE  6.  Integra!  Flux  at  10,  30,  and  100  M«y  and  Corresponding 

Characteristic  Rigidity  P  (Cont'd) 

o 


( 


Date 

J(>  10  Mev) 

J(>  30  Mev) 

J(>  100  Mev) 

P 

0 

5/  4/60 

1. 2  x 

107 

2 

(protons/cm  ) 

6  x  10* 

1. 2  x  106 

(Mv) 

127 

5/  6/60 

- 

4  x  10® 

"  c 

- 

5/13/60 

1. 5  x 

107 

4  x  10® 

4. 5  x  10 

94 

6/  1/60 

- 

4x  10® 

- 

- 

8/12/60 

- 

7 

6  x  10^ 

“  Z 

- 

9/  3/60 

9  x 

107 

3.5  x  10', 

7x  10® 

127 

9/26/60 

2  x 

109 

2.  Ox  10® 

1. 2  x  10^ 

2. 5  x  10® 
1.2x10® 

73 

11/12/60 

4  x 

109 

1.3  x  10* 

7. 2  x  10® 

124 

11/15/60 

2. 5  x 

108 

108 

114 

1 1/20/60 

1.4  x 

4. 5  x  lo' 

8  x  10 

118 

7/11/61 

1.7  x 

o  o 

0  CO  Vi 

3  x  10^ 

2. 4  x  105. 

81 

7/12/61 

5  x 

4x  10* 

3x  10® 

1  x  10? 

56 

.7/18/61 

1  X 

107 

107 

4  x  lo' 

102 

7/20/61 

1.5  x 

5  x  10® 

9x10® 

120 

9/28/61 

5  x 

107 

6  x  10 

1.1  x  10° 

121 

11/10/61 

- 

- 

- 

- 

2/  4/62 

- 

5 

5 

4 

- 

10/23/62 

6  x 

10° 

1. 2  x  10 

1  x  10 

83 

< 
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So  Ivina  Eg.  4  far  P  .  w  **  rsUtnin 

y  —  .  o/ 

o  lnlJ‘(>  3 0) / J‘ (>  100) j  M 

The  spectrum  was  fitted  at  ihe  30-  and  100-Mev  points  rather  than  the 

10  and  30  or  10  and  100  points  because  the  protons  below  30-Mev  kinetic 

energy  do  not  contribute  appreciably  to  the  dose  for  shielding  thicknesses 

2 

of  a  few  gm/cm  .  The  calculated  values  of  P  are  listed  in  Table  6  for 

o 

those  events  with  two  or  more  flux  values.  Thus,  the  unit  flux  dose  for 
th 

the  i  shielding  configuration,  d.,  is  determined  from  Fig.  15  for  the 
appropriate  P  ,  The  radiation  dose  for  each  solar-particle  event  is  then 
obtained  by  multiplying  the  unit  flux  dose  and  the  integrated  particle 
flux  for  that  event. 

The  problem  remains  to  determine  a  Pq  for  the  thirteen  remaining 

events  having  only  one  flux  value.  Since  it  is  impossible  to  calculate 

the  doses  for  these  thirteen  events,  we  must  either  associate  a  particular 

P  with  each  event,  based  on  some  relationship  of  the  observed  P  's  to 
o  r  o 

other  quantities,  or  be  content  with  arriving  at  the  most  probable  or 

mean  doses  from  the  events.  It  is  shown  in  Fig.  15  that  the  unit  flux 

dose  is  a  strong  function  of  Pq  fora  given  shielding  configuration.  This 

indicates  that  the  use  of  an  average  Pq  wili  not  necessarily  give  the  mean 

doses.  It  remains  to  be  determined  whether  the  P  for  a  given  event  can 

o 

be  related  to  some  other  quantity.  Variables  to  be  considered  include 
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the  fete!  flux  and  the  f!nrp  location  on  the  sun's  disk.  Figure  16,  a  plot 

of  the  calculated  P  ‘s  versus  the  integrated  flux  above  30  Mev  for  the  31 
o 

events  having  two  flux  points,  shows  no  apparent  correlation  between  flux 

and  Pq.  In  Figs.  17  and  18,  the  calculated  P^'s  are  plotted  as  a  function 

of  the  position  of  the  solar  flare  producing  the  event  on  the  sun's  surface. 

Again,  no  correlation  is  apparent.  From  these  figures,  P  is  apparently 

a  random  quantity.  To  check  this  conclusion,  the  number  of  events  having 

a  greater  than  a  given  rigidity  P  is  plotted  as  a  function  of  P  In  Fig.  19. 

Except  for  the  single  point  having  P^  -  195  for  the  February  23,  1956  event, 

the  points  fall  roughly  on  a  straight  line.  This  graph  supports  the  conclusion 

that  the  statistical  distribution  of  P  is  rectangular  between  45  and  150  Mv. 

o 

We  conclude,  therefore,  thcr  df/dP^,  the  probability  density  in  P  ,  is  a 
constant,  where 

for  45  <  P  <  1 50  Mv, 

(Eq.  5) 

for  all  other  P  . 

o 

(Equation  5  can  be  deduced  directly  from  Fig.  19  since  it  is  essentially  o 

plot  of  the  cumulative  probability  versus  P  . ) 

The  determination  of  the  mean  doses  is  now  straightforward.  If 

d.(P  )  is  the  unit  flux  dose  as  a  function  of  P  (Fig.  15)  for  the  i^  shielding 
i  o  o 

configuration  and  d.  is  the  mean  of  dj,  then 
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d. 

i 


.  CK3 


df 

dF 


d.tr  ;  ar 

i  o  o 


-O0 


d.(P  )  dP  . 
i  o  o 


vEq.  6) 


Thus,  we  need  only  integrate  the  curves  in  Fig.  15  to  arrive  at 
tho  various  d..  Table  7  lists  the  values  for  d.  obtained  in  this  manner. 

Table  8  (obtained  by  use  of  the  calculated  P^  when  available  and  otherwise 
by  use  of  Table  7)  lists  the  doses  by  event  for  nine  different  shielding  con¬ 
figurations. 

The  present  computer  shielding  program  can  handle  rather  compli¬ 
cated  shielding  configurations.  Figure  20  is  an  example  of  the  manner  in 
which  a  spacecraft  is  treated  as  a  shielding  configuration.  The  spacecraft 
is  divided  into  an  appropriate  number  of  sectors  and  the  number  of 
grams  per  square  centimeter  is  determined  for  each  kind  of  material  con¬ 
tained  in  that  sector.  This  constitutes  a  shielding  configuration  designated 
by  d,.  The  existing  program  can  handle  up  to  350  sectors  and  up  to  as 
many  as  twelve  layers  of  different  materials. 
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TABLE  7.  Mean  Unit  Mux  uose  a.  for  Nine  oiucrlJIny  Configuration; 


Shielding 

Configuration* 

d.  ** 

10/4 

(rad/proton  >  30  Mev) 

9.  89  x  10"9 

6/4 

1.67  xl0“8 

4/4 

2.35  x  10"8 

2/4 

3.54  x  10“8 

1/4 

4.32  xl0”8 

1/0 

4.  75  x  10'7 

2/0 

2.  08  x  Iff7 

5/0 

5.99  x  10“8 

10/0 

2/16  x  10"8 

*  see  Table  8 

*  see  text 
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11/12/60  484  263  100  43  75  62  46  34  20.8 

11/15/60  288  151  53  20.5  39.6  31.7  23.0  16.6  10.1 

11/20/60  17.3  9.4  3.4  1.3  5  2.52  2.07  1.49  1.08  0.68 


TABLE  8.  Doses  by  Event  for  Nine  Shielding  Configurations  (Cont'd) 
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GROSS  RADIATION  DOSE  PREDICTION 

Cumulative  orobaoility-dose  curves  for  cycle  19  are  constructed 
from  Table  3  as  follows.  The  mission,  of  duration  T  days,  is  considered  to 
start  on  each  consecutive  day  beginning  January  1,  1956  and  in  each  case 
the  last  minion  end*  on  Oecember  31,  1962.  For  each  starting  date,  the 
dose  which  would  lwve  been  received  is  recorded,,  After  all  the  doses  have 
been  calculated,  they  are  arranged  in  descending  size.  The  number  of 
mission  intervals  which  sampled  a  dose  larger  than  a  given  dose  D  may  now 
be  determined  by  numbering  the  descending  doses  consecutively.  The 
probability  of  receiving  a  dose  greater  than  D  is  then  the  ratio  of  this 
number  to  the  total  number  of  missions  (25 57  -  T).  (There  are  2557  days 
from  January  1,  1956  through  December  31,  1962.) 

Representative  curves  for  mission  durations  of  thirty  and  sixty  days 

2 

with  shielding  thicknesses  of  1,  5,  and  10  gm/cm  of  aluminum  are  shown 
in  Figs.  21  and  22,  A  program  which  carries  oui  this  portion  of  the 
calculation  has  been  written  for  the  IBM  7090.  This  program  supplements 
the  existing  computer  programs  for  calculating  the  doses  and  allows  rapid 
construction  of  a  probability-dose  curve  for  any  desired  mission  duration 
and  shielding  configuration. 

In  a  similar  fashion,  cumulative  probability-particle  flux  curves 

may  be  derived.  Figure  23  gives  the  probability  of  receiving  an  integrated 
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particle  flux  (E  >30  Mev)  .greater  than  a  given  integrated  flux  J  for 
exposure  times  of  3,  6,  14,  and  24  months. 

The  previous  section  of  this  report  dealing  with  solar-particle  event 
predictions  indicates  that  the  probability  of  encountering  solar-flare-particle 
events  can  be  substantially  reduced  by  postponing  scheduled  missions  or 
aborting  missions  in  progress  if  certain  solar  conditions  occur.  The  extent 
of  the  reduction  and  the  approach  to  be  employed  will  expend  upon  the 
mission  type  and  mission  duration. 

Establishment  of  solar  conditions  to  be  used  and  calculations  of 
the  reduced  probabilities  for  various  missions  is  now  under  tray  in  the 
Space  Physics  Group  of  The  Boeing  Company. 

Some  information  concerning  solar-particle  events  is  summarized  in 
Fig.  24,  which  is  a  plot  of  the  smoothed  sunspot  numbers  for  solar  cycles 
18  and  19  with  the  maximum  and  minimum  values  predicted  for  solar  cycle 
20  (Refs.  21  end  23).  The  number  of  class  3  and  3+  flares  per  year  and 
the  number  of  particle  events  per  year  are  also  shown  in  this  figure.  The 
number  of  class  3  and  3+  flares  for  cycle  20  are  predicted  on  the  basis  of 
correlations  derived  from  cycles  18  and  19  (Ref.  24). 

Examination  of  Fig.  24  indicates  that  the  use  of  probability  curves 
obtained  for  solar  cycle  19  would  apply  to  solar  cycle  20  if  solar  activity 
reaches  the  maximum  predicted.  The  dose  probabilities  obtained  from 
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f 

such  curves,  than,  would  certainly  bo  an  upper  limit  for  doses  to  He 
Incurred  during  corresponding  missions  in  cycle  20.  Better  predictions 
concerning  the  sbapn,  the  magnitude  of  the  peak  ond  time  of  peak  of  the 
smoothed  sunspot  curve,  with  corresponding  changes  in  the  predicted  dose- 
probability  curves,  can  be  made  after  the  beginning  of  solar  cycle  20 
(Rot.  21). 
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